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Abstract 

We consider operators of Kramers-Fokker-Planck type in the semi-classical 
limit such that the exponent of the associated Maxwellian is a Morse func- 
tion with two local minima and a saddle point. Under suitable additional 
assumptions we establish the complete asymptotics of the exponentially small 
splitting between the first two eigenvalues. 

Resume 

On considere des operateurs du type de Kramers-Fokker-Planck dans la 
limite semi-classique tels que I'exposant du maxwellien associe soit une fonc- 
tion de Morse avec deux minima et un point selle. Sous des hypotheses 
supplementaires convenables on etablit un developpement asymptotique com- 
plet de I'ecart exponentiellement petit entre les deux premieres valeurs pro- 
pres. 
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1 Introduction 



This paper is a natural continuation of the work [15], investigating the low lying 
eigenvalues of the Kramers-Fokker-Planck operator 

P = yhd,-V\x)-hdy + ^{-hdy + y)-{hdy + y), x,i/eR", (1.1) 

where 7 > 0. Physically the semiclassical limit h ^ Q corresponds to the low 
temperature limit. As explained in [15], the original motivation for that work was 
to give more explicit versions of some results in [H] and later in [9] , giving estimates 
on the time of return to equilibrium, or more or less equivalently, on the gap between 
the first eigenvalue (when the potential V tends to plus infinity sufficiently fast at 
infinity) and the second eigenvalue. See also [25] for further developments in that 
direction. The methods of those works as well as the one of Eckmann and Hairer 
[H] are inspired by those of hypoellipicity for Hormander type operators. We will 
not repeat here all the motivations of [15], coming also from the works [H [17] and 
others, including recent developments in pseudospectral theory. 

In [15] it was assumed that V in (11.11) is a smooth real Morse function with finitely 
many critical points, t/j, j = 1, .., A^, that d"'V{x) is bounded for all multiindices a 
of length |a| > 2 and that |Vy(x)| > 1/C for \x\ > C for some C > 0. Under these 
assumptions it was shown that the eigenvalues in any disc D{0,Ch) are of the form 

X,4h) ~ /i(/i,- + + /i'/^-^/i,,fc,2 + ...), h ^ (1.2) 

where the index j labels the critical points and the leading coefficients fij^k can be 
given explicitly in terms of the Hessian of V at the corresponding critical point. The 
values fij^k are confined to a sector {|arg(2 — /ij,o)| < for some 9j E [0, |[, and 
Aj^o is the eigenvalue with the smallest real part of all the Xj^k- (This comes formally 
from a harmonic oscillator approximation.) Further /Xj ^ 7^ /ij for A; 7^ and the 
asymptotic expansion (11. 2p for Xj^ contains only integer powers of h. See |15j for 
more details. Using this result, as well as control over the resolvent along suitable 
contours in the right half plane, the authors were able to give asymptotic expansions 
for large times of exp {—tP/h), that emphasize the role of the eigenvalues close to 
given in (11.21) . Indeed, there are no others in a certain parabolic neighborhood of 
the imaginary axis. 

Moreover, we have fij^ > with equality precisely when Uj is a local minimum 
of the potential, and in the case of such a minimum it follows from the above results 
that Aj^o is actually exponentially small. 

In this paper we address the question of determining more precisely the size of 
these exponentially small eigenvalues. In the case when V has precisely one local 
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minimum, say Ui, and tends to +00, when x 00, we know that the corresponding 
eigenvalue Ai^o is equal to zero (with the Maxwellian exp (— (?/^/2 + V{x))/h) as the 
corresponding eigenfunction) and that this eigenvalue is separated from the other 
ones by a gap of size h. This means that we have return to equilibrium with a 
speed that is roughly 1. The situation becomes more complicated when there is 
more than one local minimum. We are then in the presence of a tunneling problem 
which is much more complicated than the corresponding ones for the semiclassical 
Schrodinger operators since our operators are non-elliptic. In principle one should be 
able to follow the general approach of earlier works in the Schrodinger case as [TO] . 
However it seems that one necessarily runs into a tunneling problem where the wave 
functions have to be studied also in a neighborhood of some intermediate saddle 
points of V, and as known from [llj that can indeed be done in the Schrodinger 
case with techniques that are very useful in a variety of problems. To carry out such 
an approach in the case of Kramers- Fokker- Planck would require one to accumulate 
the difficulties of non-resonant wells with the ones coming from the lack of ellipticity. 
This seems to lead to considerations of degenerate non- symmetric Finsler distances 
(see for instance [Tl [TU]). 

For the Witten Laplacian (see [12]) we are also in the presence of a tunneling 
problem with intermediate non-resonant wells and in that case one could avoid 
the detailed study near the non-resonant wells by studying directly the Witten 
complex as a tunneling problem between critical points of neighboring indices. More 
recently M. Klein, B. Helffer, and F. Nier [8] have used that approach to study the 
exponentially small non- vanishing eigenvalues of the Witten Laplacian. Also in [H], 
explicit estimates relating such small eigenvalues for the Witten Laplacian and the 
Kramers-Fokker-Planck operator were established. 

This relation with the Witten complex was strengthened further in the works 
of J. Tailleur, S. Tanase-Nicola, J. Kurchan [22] and J. M. Bismut [2J, who showed 
using respectively the languages of supersymmetry and differential forms, that the 
Kramers-Fokker-Planck operator can be viewed as a Witten Laplacian in degree 
associated to a certain non-semidefinite scalar product in the spaces of differential 
forms. See also [18] for a quick introduction to the differential form version of Bismut 
and [3]. 

In the present paper, we use this supersymmetric approach. Our main result, 
valid also for a class of more general operators, is that if a certain weight function 
(which in the KFP-case is the function + V{x)) has precisely two local minima 
U± and an intermediate saddle point Uq then we can get a complete asymptotic 
expansion for the second eigenvalue of the corresponding Witten Laplacian (reducing 
to the KFP operator in the special case). The logarithm of this eigenvalue is equal 
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to —2h^^{mm{(f){Uo) — (f>{Ui) , (f){Uo) — 0(f/_i)) + o(l)), but actually we do have a 
complete asymptotic expansion. See Theorem 111.11 for a complete statement. It 
seems clear that this result can be somewhat generalized but a more complete result 
might require exponential estimates and asymptotics for eigenf unctions also far from 
the critical points. In our present approach we are able to get such information 
for the eigenfunctions of the degree operator in the basin of attraction of each 
minimum and for the degree 1 Laplacian in a small neighborhood of the saddle 
point. 

In most of the paper we work with a scalar real second order non-elliptic operator, 
which is also non-self adjoint, and we were led to reconsider some steps in [15]. The 
plan of the paper is the following: 

In Section [2] we do some very simple and elementary exponential estimates mainly 
designed to get the appropriate control near infinity. 

In Section [3] we establish the m-accretivity for our operators so that the step from 
a priori estimates to spectral information becomes possible. 

In Section m we study certain auxiliary weights, somehow related to escape functions 
in resonance theory (see [13] and a large number of more recent works) in connection 
with some dynamical conditions. 

In Section O we use those weights together with a machinery of Fourier integral 
operators with complex phase in order to get phase space a priori estimates away 
from the critical points. This section is perhaps technically the most complicated 
one, but the underlying ideas are now quite standard. Alternative methods are 
certainly possible and we might return to this step in future works. This section 
and the subsequent one are quite technical and should not be studied in detail in 
the first reading. 

In Section [6] we study the conjugation of our original operator under the Fourier 
integral operators of the preceding section and in Section [7] we finally obtain the a 
priori estimates that we need. 

In Section [S] it is now quite easy to get detailed asymptotic results by adapting the 
methods of [15j. 

In Section [9] we show that the eigenfunction associated to Xj^ for each (nondegen- 
erate) critical point Uj has the form a{x; h)e~'^^^^^^ in a neighborhood of that point, 
where a{x; h) has an asymptotic expansion in integer powers of h with coefficients in 
C°°(neigh [Uj)) and (p = (pj is a smooth function of the order of magnitude \x — Uj\^. 
It would be very interesting to extend such descriptions further away "beyond caus- 
tics" . 

In Section [TU] we review the supersymmetric approach of [22], l2j (see also [IB]) and 
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establish various interesting links between the dynamical conditions of Section [Hand 
old results for non-self adjoint operators with double characteristics ([23]. This sheds 
additional light on some related computations in [T^ . 

Finally in Section [TTl we can put the various results together and establish the 
precise exponential asymptotics of the spectral gap between the first and the second 
eigenvalue (both real and the first one being zero.) 

We expect that the spectral results of the present paper will give rise to precise 
asymptotics for the associated heat-evolution problem in the limit of large times and 
we plan to treat that problem in a separate paper. 



Acknowledgements. We are grateful to B. Helffer who pointed out the work 
[22j to us, and for pointing out that an earlier version of our condition (11.2) was 
not directly applicable to the Kramers-Fokker-Planck operator. The research of 
the second author is supported in part by the National Science Foundation under 
grant DMS-0304970 and by an Alfred P. Sloan Research Fellowship. He is happy to 
acknowledge the hospitality of Ecole Polytechnique and Universite de Reims, where 
part of this work was done. 



2 A priori estimates 

In this section we establish some simple a priori estimates which will be important 
in Section E] and at other places. They illustrate the technique of gaining ellipticity 
by means of exponential weights that we shall later employ also in a micro-local 
setting. 

Let M denote either the space R*^, or a smooth compact n-dimensional manifold 
equipped with a strictly positive smooth density of integration dx. On M we consider 
a second order differential operator 

n 1 " 

P = ^ hD^^ o hj^k{x) o hD^^ + - ^{cj{x)hd^^ + hd^^ o Cj{x)) + Pq{x) (2.1) 

j,k=l j=l 

I OXj 

where the coefficients bj^j^, Cj, po are assumed to be smooth and real, with bj^k = bkj- 
In the manifold case, we use local coordinates such that dx = dxi...dxn- To P we 
associate the symbol in the semiclassical sense, 

p{x,i) = P2{x,i) + ipi{x,C} +po{x), (2.2) 
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P2{x,0 = ^hk{x)^jik, Pi{x,0 = 'Ycj{x)^j, (2.3) 

j,k=l j=l 

SO that pj{x,^) is a real- valued polynomial in ^, positively homogeneous of degree j. 
(It is well-defined on T*M and coincides with the Weyl symbol mod 0{h?) locally 
uniformly.) We assume that 

P2(a:,O>0, Po{x)>Q. (2.4) 

In the case M = R", we impose the following growth conditions at infinity: 

d:b^4x) = 0(1), |«| > 0, (2.5) 

d^cjix) = 0{1), \a\ > 1, (2.6) 

d^poix) = 0(1), |a| > 2. (2.7) 

When discussing P in the operator theoretical sense we will assume that it is the 
closure of P : S{M) S{M) as an unbounded operator in L'^{M). (When M is 
compact, we identify the Schwartz space S{M) with C°°{M).) Let V{P) C L'^{M) 
be the domain. 

Lemma 2.1 We have 

j Po{x)\u{x)\^dx + jj2Kkix)ihD^^u)(hD~il)dx = Re{Pu\u), (2.8) 



for allu e S{M). 



Proof Immediate by integration by parts. □ 

In the manifold case, we view B{x) = {bj^k{x)) as a positive semi-definite matrix 
T*M — ^ TxM and if we choose some smooth Riemannian metric on M, we can view 
B{x) as a map T*M — > T*M and define B{x)2 similarly. (12.81) then becomes 

\\po{xf^uf + \\B{x)hDuf = Re (Pn|n), (2.9) 

implying 

\\pIu\\ + \\BhDu\\ < Coi\\Pu\\ + \\u\\). (2.10) 
In particular, u G T>{P) =^ ||Po^|| + ||i?2/;,Z}n|| < oo. 
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Using the anti-selfadjoint part iPi we shall obtain a similar estimate where the 
averages of po along the trajectories of 

n 

iy{x,d^) = ^Cj{x)d^. (2.11) 

will play a role. 

In general, if iplx) is a smooth real-valued function, the operator 

:=e^/'^oPoe-'^/'^ (2.12) 

is of the same form as (12.1 1) with new coefficients Cj, Pq and the new symbol 

= P2ix,^ + iij'ix)) + ipi{x,^ + i^lj'{x))+po{x) (2.13) 
= P2(a;,0 +«(Pi(a;,0 + d^p2ix,ip'{x)) ■ ^) 

+Po{x) - pi{x,tjj'{x)) -p2{x,^'{x)) 
=■ P2ix,^) +ipi,^{x,^) +po,^{x) 

In this section we choose very small and treat p2{x,ilj'{x)) as a perturbation. 
Notice that 

p^ix,^'{x)) = ly{x,^^)^p. (2.14) 

Let f{t) e C°°([0, oo[; [0, 3/2]) be an increasing function with f{t) = t on [0, 1], 
f{t) = 3/2 on [2, oo[, f{t) < t. Put f,{t) = ef{t/e) and consider for Tq > fixed. 



j k{^)feOpQoexp{ti^)dt, (2.15) 



where 

Then 
where 

{feOpo)To = rp 

J-O J -To/2 

is the time To average of o p^ along the integral curves of u. Clearly, 

Tr\ 3 3Xn 

\ipe\< —-^ = —e. 2.18 

' - 4 2 8 ^ ' 



0, \t\ > 1/2, 
k{t) = <( t + |, -| < t < 0, 
-k{-t), < t < i. 



^ii^e) = feOPO- {feOPo)To, (2-16) 

1 rTo/2 

{fe o Po)to = ttt / /e o Po o exp (tz/)dt (2.17) 
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From (12.61) it is easy to see that := exp (th'){x) is well-defined for all t G R, 

X G M and that 

< C(a)e^ol"ll*l, a G N", |a| > 1. (2.19) 
In particular, |c}"$f(x)| < Ca,To, \t\ ^ ^o- On the other hand, since po > 0, = 

C(l) (by ([22D), we know that \p'o\ < 0{l)p^ and this quantity is 0{e^/^) in the 
region where < po < 2e. It follows that 9°(/e o p^) = C(e^~l'^l/^), a G N" and 
together with (12.191) . we get 

= C(e'-l"l/2), Va G N", (2.20) 

for every fixed Tq. 

Then from (|2J[3D . fl27Hl) . fl236D . (K20\\ and the uniform boundedness of the 6^-^, 
we get 

Po,5i/;. = Poix) - 6f,opQ + 6{f,opQ)j.^ - 6'^p2ix,tlj'^) (2.21) 
= Po(x) - 6f,opQ + opo)ro - C'('^^e), 

uniformly on R" for < 6 < 1. Using the properties of /, we notice that po{x) — 
SfeOPo > (1 - S)po, so 

Pom. > (1 - ^)Po + HfeOPoWo - 0{6h). (2.22) 

Since the coefficients of P^^^ grow at most polynomially. Lemma 12.11 can be 
applied and gives 

j (PomA^) -Re2;)|Mp(ix < Re{{Ps^, - z)u\u), z e C, Rez < (2.23) 

where Z ^ 1 is independent of e. (For Lemma [2TT] we do not need that Po > 0.) Let 
H > and rewrite (12.231) as 

J max{po^^, ij,)\u\'^dx < Re (((P^ - z) + {jj, - po^^)+)u\u) 
= Re (max(po,v>,yu)~^(Pv. - z + {n - po,rp)+)u\m8ix{po^^, fi)^u). 

Here we write for 6ipe and po,^ = po,ip — Re z. 
Then using Cauchy-Schwarz, 

II /~ \i II II /~ \ — i/n \ II II POl/')+ /~ \i II 

II max(po,v,/^)^M|| < ||max(po,v,W H-^^/' " -^J^ll + II ^ — max(po,v, Ai)^M||, 
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leading to 

II max(po,v,/^)^M|| < II max(po,^,/i)"5(p^_2;)M|| + ||l{pQ_^,<^} max(po,,/,,/^)^M||. (2.24) 

Notice that l{po ,^</,} max(po,v, A*)^ = ^Po,i,<^^}^^^ ■ 
We have 



> Po{x)-6f,op, + 6{f,op)T^-0{6'e) 



Choose Z = 6 ^, so that 

Po,i^ > Po{x) - Sf.opQ + S{f,op)T^ - 0{5'^e). 

Choose 6 > small enough (but independent of e) so that 0{6^e) < 5e/Co, where 
we shall fix Cq sufficiently large. Then 

PO,i^ > (1 - S)Po{x) + S{{feOpo)To - 7^). 

When (/e opo)to > f , we get 

Po,v > (1 - S)Po{x) + —■ 

When (/, o po)ro < f we have po,^ > (1 - 5)Po - ^■ 
Choose 

(5e 

Then 

(1 - 6)po + li < Po,i> + 2/i < 3 max(po,v, A^)- (2.25) 

Moreover, if pq,^ < fi, then 6{f,opo)To - ^ < J, so (/e opo)ro < ^- From ([221 
and (12.251) we then infer that 

\m-S)po + ^t^u\\< (2.26) 

I Sc 

3||((1 - 6)po + -^yHPi' - z)u\\ + VsJ 7rlkll{(/<opo>Tn<^}' 

Oq V '-'0 " 



for Re z < -|. 
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Here we can take e = Mh with M ^ 1. Then by (12.181) we have 

\i,/h\ = \5i,Jh\<C{M,T,)5, 
so there is a constant C independent of h (but depending on M, 5) such that 

From the discussion above, in particular (I2.26p . we get 

Proposition 2.2 Let P he of the form ( 12. ip . where bj^k,Cj,po are smooth and real 
and satisfy (I2.2p - (12.71) . Define {fe ° Po)to as (I2.17P with f^ defined after (I2.14p . 
Then for every C > 0, there exists C > such that 



50 + h)-^u\\ < C(||(po + hyHP - z)u\\ + /^^lh^ll{(/g^opo>To<cM)' (2.27) 
for M G iS, Re 2 < Ch. 

Notice that (12. 8p imphes that 

\\BhDuf < \\{po + h)-^{P~z)u\\\\{po + h)^u\\ + C\\h'^uf. (2.28) 

3 From injectivity to the resolvent 

Let P = P2 + iPi + Pq with symbol p = p2 + ipi + Po be as in section [21 so that we 
have 

Po{x)\u{x)\'^dx + {P2u\u) = Re{Pu\u), u e S{M), 



leading to 

|Re;z|||Mf + \\pluf + {P2u\u) = Re{{P-z)u\u) < \\{P - z)u\\\\u\\, Rez < 0. 
We get 

||(|Re2;| +Po)-^uf < \\{\Rez\ + po^kP - z)u\\\\{\Re z\ +Po)-^u\\, 

\\{\Rez\+po)-^u\\ < \\i\Rez\+po)-HP-z)ul 
iRe^l^lKlRe^l +Po)^u\\ < \\{P - z)u\\, 
|Re2|||n|| < \\{P-z)u\\ 
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From this we get 

{P2U\U) < \\{P-Z)U\\\\U\\ < -J—\\{P~z)uf, 

|Ke z\ 

and putting some of the estimates together, 

\Rez\'^\\uf + \Rez\\\pluf + \Rez\{P2u\u)) < 2\\{P~z)uf. (3.1) 

By P we also denote the graph closure of P : S ^ S. From the estimates above 
we see that the range 71{P — z) is closed in when Re 2; < 0. 

Proposition 3.1 

n{P-z) = L^, Rez<0. 

Proof It suffices to prove that TZ{P — z) = for some z with Re z < 0, because the 
a priori estimate then implies that ||(P — z)^^\\ < |Re2;|~^ and this fact extends by 
standard arguments to the whole left half plane. For the same z it suffices to show 
that if M G and (P* — 'z)u = in the sense of distributions, then m = 0. Now 
the formal adjoint P* = P2 — iPi + Pq has the same properties as P, so in order to 
simplify the notations, we may just as well prove the corresponding fact for P — z in 
stead of P* — z: There exists a z with Rez < such that if u G and [P — z)u = 
in the sense of distributions, then u = 0. 

When M = R", let Opf^{q) denote the Weyl quantization of q{x, h^) and put 

Ae = Op,((e(x,0)""^) 

where iV > 2 is fixed and e > is small and fixed. (When M is a compact manifold, 
choose a Riemannian metric and put = (1 — e/i^A)"^/^.) Consider the equation 

(P — z)u = V, u,v E 

Then A,(P - z)u = A,w, so 

(P - z)Ku = A,v + [P, A,]u. (3.2) 

1 1 

Since > 2, we can find a sequence uj G S such that A^uj — > A^u, PqA^uj — > p^A^u, 
Ae(P - z)uj A^v, [P,A^]uj [P,A^]u in L^, {P2A^Uj\A^Uj) (P2A,u|A,u). This 
means that fl3.ll) is applicable to fl3.2p . and we get 

\Rez\^\\A,uf + \Rez\\\plA,uf+\Rez\{P2A,u\A,u) < 4(||A,t;f + || [P, A,]nf ). (3.3) 
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Consider [P,A^]u = {[P, A^]A-^)A^u. We study the operator [P, A.jA^^ and 
assume that M — R" in order to fix the ideas. If p = (x, we have 

a,"(p)-^ = o(i)(p)-^-i"i, 

so 

a;(6p)-^ = ompr'^i^p = 0(i)(6p)-^(p)-i"i, 

uniformly with respect to e. Prom this we deduce that the symbol of [P, Ag] is 
Ag} + Oo(^^Ag), using the same letters for operators and their symbols (except 
for P, Pj where we already introduced a distinction by using lower case letters for the 
symbols), and using the notation Oo{h'^A^) for a symbol q satisfying d^q = Oa{h?A^), 
uniformly in e. Here 

{p,AJ = Oo(A,(p)(p)), 

so the symbol of [A^, P]A^^ is 

Recall that p{x, — Po(3^) + +^2(2^, 0- We get 

{pi,AJ _ Ax). (^^h-n (^) 
= g + c»o(i), Q = • = a(^) • vMp)- 

From these computations, we retain that 

\\[P,A,]A-;^w\\<0{h)\\w\\+0{h)\\Qw\\. (3.4) 

The symbol Q is real and d^Q = Co((0^"'"'/(p))> so 

||g^/;||2 = (Q2^|w), 

where Q oQ has the symbol + Co(^^)- 
Since p2 > 0, dlp2 — C(^^) we know that 

\d,P2\<O{imP2{x,0K 
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and we conclude that on the symbol level 

Hence by the semi-classical Fefferman-Phong inequality for operators with symbols 
in the Hormander class 5*^ q (see for example Subsection 7.2 in [15] for a short review 
of the Weyl-Hormander calculus): 

WQwf <C{P2w\w) + 0{h'')\\w\\\ 

Using this in ( 13. 4p . we get 

||[P,Ae]A7'wf < Ch\{P2w\w) + \\w\\''). (3.5) 

we use this in (13.31) with [P^k^u = [P^K^K'^^keU and get 

\Re z\\\pIk,u\\'^ + {\Re - Ch^)\\Ku\\'^ + {\Re z\ - Ch^){P2Ku\Ku) (3.6) 

<4||A,(P-2)m||1 

So if Re z < —VCh and (P — z)u = 0, u E L^, we have u = 0. □ 

Corollary 3.2 The maximal closed extension P^aax of P (with domain given by 
{m G L^; Pu G L^} coincides with the graph closure (the minimal closed extension), 
already introduced. 

Proof Let z be fixed with Re^; < —\/Ch. Let u,v G with (Pmax "~ z)u = v. 
Denote by Pmin the graph closure. Since 7^(Pmin — z) = L^, there exists u G ©(Pmin) 
such that (Pmin — z)u = V. Hence {P — z){u — u) = in the sense of distributions. 
We saw in the proof of the proposition that this implies that u — u = 0. □ 

This result can be extended to the following auxiliary problem: Let P_ : 
L^, P+ : —>■ be bounded operators and assume that 

Viz) = ^ ) : V{P) X ^ X C^, 

is injective with 

\\u\\ + \u^\ < Cih){\\v\\ + \v+\), z G n{h), (3.7) 



whenever 



n^) ( „n = l!" ) • (3-8) 
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Here we assume that fl{h) is open and connected, intersecting the resolvent set p(-P) 
of P. 

From (13. 7p . we see that V(z) is injective and has closed range for every z G fl{h). 
When z is also in the resolvent set of P, we can rewrite (13.81) as 



1 (P - z)-^R_ \ f u \ _ f {P- z)-^v 



(3.9) 



and the matrix appearing here is injective. On the other hand it is a finite rank 
perturbation of the identity in x and the injectivity implies the bijectivity. 
It follows that V{z) is bijective for z G Q{h) fl p{P). Combining this with (13.71) and 
recalling that Q is connected we see by a standard argument that V{z) is bijective 
for aWzen and that \\V{z)-^\\ < C{h). 

4 Geometric preparations 

In this section we construct a special bounded weight that we shall implement in 
Section [5] with the help of Fourier integral operators with complex phase. 

Let P be the operator introduced in the beginning of Section [2l satisfying (12. ip - 
flO) . The symbol 

p{x,0=Po{x) + ^^^ (4.1) 

is non-negative and satisfies 

{d,A0d^rP = O{l), \a\>2, (4.2) 
locally uniformly with respect to x. It follows that 

{d,A0d^rP = Oip-^), \a\ = l. (4.3) 
We now introduce a critical set associated to p. 
Hypothesis 4.1 Assume 

The set {x G M; po{x) = 0, h'{x,dx) = 0} is finite = {xi, ...,xn}- (4.4) 
Let pj = {xj, 0) and put 

C = {pu-,Pn}. (4.5) 
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Notice that Pi,Po-,P2,P vanish to second order at each pj. Our weight will be of the 
form 



- j kijr)Pe o exp itHp,)dt, < e < 1, (4.6) 

where Pe will be specified below and k is the same function as in Section [21 Notice 
that 

Hp,A = {Pe)To -Pe, (4.7) 

where we now write in general, 



1 /•To/2 

{(1)to = 7^ g o exp {tHpJdt. 

J-O J -To/2 



Let g{t) G C°°([0, +oo[; [0, 1]) be a smooth decreasing function with 

, , fl, < t < 1, ,^ 

Notice that g^^\t) = 0{{t)-^-^). In a domain < |p - Pj| < 1/C, C > 1, we put 

P;(p)=^7(^^^)p, (4.9) 
e 

so that P'e(p) = p{p) for |p — Pj| < \fe. It is easy to show the symbol estimates 

d'-p. = g( (^^|^_'^^.p)H/. )> < |p-p,| < l/C. (4.10) 

Away from any fixed neighborhood of C, we simply get 

d"p, = 0{e). 

Let Xj ^ (^^(R^") be equal to 1 near pj and have its support close to C. To define 
further out from C, we put 

Pe(p) = 5^X,(p)^7(^^)p + e(l-5^X.(p))p, (4.11) 
1 

so that (14.101) remains valid in a neighborhood of C, while 

a^afp, = O{e{0-^^\) = 0(6(0-1^1) (4.12) 
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outside any such region where |x| is bounded. We also have 

d:dl{i-d^p;) = 0{e{i)-'-\^\) (4.13) 
in the same region. In fact, fl4.13p follows if we write for |,^| ^ 1: 



From (14.111) we see that Peip) = when dist (p, C) > 1/C and in particular for 
p = (x, ^) with — > 1/C, Vj. Further away from {xi, ...,xn} we want to make 
Pe{x,^) independent of ^, so we replace p = po{x) +^2(2^, 0/(0^ (14.111) by 

Pnew(a;, = Po{x) + ^P2(x, 0, (4-14) 

where x ^ C^{M; [0, 1]) is equal to 1 near {xi, ...,X]y}. 

Even further out (in the case when M = R") we want to avoid problems caused 
by po being large, so when 2> 1 we want to replace Pe = epo there by fe{po), 
where /^(t) = ef{t/e) is the function introduced in Section [2l Thus with a new 
cutoff function Xnew G C*^(M; [0,1]) being equal to 1 in a large neighborhood of 
supp X, we get the final choice of p^: 

Pe{p) = Xnew{x)(^Xj{p)9i- —)p + e{l - ^Xj{p))Pne^) (4-15) 

1 ^ 

+ (1 -Xnew(x))/e(po)- 

Notice that by construction 

Pe<P<P0+P2- (4.16) 

We also get 

Proposition 4.2 We have (14.101) near C and (14.121) . (I4.13p in any closed subset 
of T*M disjoint from C and supp (1 — Xnow) x R-*^; where Xnew is the wider cutoff 
in (14.151) . Over supp (1 — Xnow) (where p^ only depends on x) we have d'^p^ = 
C(ei-l-l/2). 

From the definition (14. 6p . we see that ipe satisfies the same estimates as Pe and 
will only depend on x for x outside a neighborhood 7ix{C) (that we can choose as 
small as we like) and that the region where we only have d^ipf,{x) = (^(e^"'"'/^) can 
be any neighborhood of infinity. This follows from the fact that our various symbol 
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estimates (as well as the ^-independence) are conserved by the flow of iJp^, as for 
(I4.13p . we here also use that H^^ commutes with dilations in ^. 
It also follows from the construction that 

Pe{.p) = p(p), when \p - pj\ < y/e, (4.17) 
p,{p) - ^-^p, when <\p- Pj\ <1/C, (4.18) 

\p-pjr 

Pe{p) ^ epnew, whcu dist (p, C) > 1/C, Ip^I < C, (4.19) 

Pe(p) = e/e(Po), when |p^| > C. (4.20) 
We introduce the following dynamical conditions where To > is fixed: 

Hypothesis 4.3 

Near each pj we have (p)to > p^lp ~ PjP; (4-21) 

In any set |x| < C, dist (p,C) > (4.22) 
we have (p)to(p) > 5(C) > 0. 

When M = Ti^ we also need a modified dynamical assumption 

V neighborhood U of vr^C, and V a; G R" \U, 3C > 0, (4.23) 
meas({t G [-y,y]; po(exp tz/(x)) > ^}) > ^. 

When M is compact, we just assume f l4.2ip . f l4.22p . where it is understood that the 
estimate in f l4.22p should hold for all p G T*M of distance >l/C from C. 

Notice that in the region dist (x,7r^(C)) > l/C, \x\ < C, fl4.22p is equivalent (up 
to a change of C{C)) to 

(Po)to,.(x) > CiC) > 0, (4.24) 

C (C j 

where 

-.To/2 



(?)to,i' = 7^ / g o exp {tu)dt. 

Jo J -To/2 
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This follows from the fact that o exptHp^ = exp (tu) o tc^, where vr2.((x, ^)) := x 
and the fact that = 0} is invariant under the Hp^ flow. 

Actually, if we only assume (14.211) . then we see that {p2)To{xj,0 ^ and if we 
also assume (14.241) . we recover (I4.22p . first over a neighborhood of each xj and then 
again by (14.241) over any bounded set in M. 

Assuming the dynamical conditions flOTD . KTI^ . ( KMi . (K^ ). we see that 

(p,)to - dist (p, Cf, dist (p, C) < v^, (4.25) 

(p,)to -e, dist(p,C) > v^. (4.26) 

From the estimates on ip^ in the various regions that we mentioned after Propo- 
sition 14.21 we shall often only retain that 

d^d^Mx,0 = C(e^-l°+'5|/2(0~'^'), (4.27) 
and we write this for short as 

Similarly, we have in view of (14.13p . and ( 14.10p also valid for ip^, that 

e -95^, = 0(6(0-'). (4.28) 

From [21] let us recall that if / is a function locally defined on complexified 
phase space, then at every point where df = 0, we have 

Tj rrRecr Tjlma fj j tt rrRecr Trim a (a r)r\\ 

- ^Ref - ^Im/' ^if " ^ f " ^ ~lm f " ^Kefi (4.29) 

where 



H Y^A_^A) 

^^d^jdxj dxjdij 

is the complex (1,0) Hamilton field of / with respect to the complex symplectic 
form cr = ^ d^j A dxj and Hj = Hf + Hf is the real vector field which acts as 
Hf on holomorphic functions. H^^J denotes the real Hamilton field of Im / with 
respect to the real symplectic form Im a, and similarly for the other Hamilton fields 
appearing in (I4.29p . As usual, J denotes the action on tangent vectors induced by 
multiplication by i. (When M is compact we may assume without loss of generality 
that M is real- analytic.) 

Assume first that M = R" and put Aq = r*R" = R^", 

A, = {p + tsH^^ (p); pe Ao}, < s < 1. (4.30) 
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If we extend V'e to be a function on the complexified phase space C^", by setting 

Mp)^M^p), (4.31) 
then we have the equivalent description of as 

A, = exp(si/J-'^)(Ao). (4.32) 

It follows that As is an I-Lagrangian manifold, i.e. a manifold which is Lagrangian 
for Im a. 

Now we would like to parametrize A^ by means of a canonical transformation 
k{s) : Aq — As for the symplectic form Re a (implying that A^ like Aq is an IR- 
manifold, i.e. a manifold which is Lagrangian for Ima and symplectic for the re- 
striction of Recr). For < s -C 1, let ipe,s = ^{s, p) be the unique function which is 
affine in Im p and satisfies 



d 



dp'ipe,s = on As, i>e,s = A ou As. (4.33) 

/ ^ d „ 



Recalling that V'e(p) = '^e(Rep), we get with the notation {ip^)'^ = -r^^ipe, {i'eyx ~ 



SRe X 



ipe to indicate derivatives in the real directions. 



^e,s{^^ = i^e{x, ^) + iax{Re (x, 0) -{Imx- s{iIjJ^) + ia^{Re (x, 0) • (Im^ + s(^/',)^), 

(4.34) 

where 



l-ts%X^e{x,0)) (^j) = (l^j?) (4.35) 

(4.36) 



and 



is the fundamental matrix of It also follows from the construction that 

a, = a(e^), a^ = a(e^(0-'). (4.37) 



Since ip^^s = V'e on A^, we know that — H^'^ is tangent to As for every 

s. We can therefore define «;(«) : Aq — > As, by 

^A'){p) = H'ZiM'W)^ P e ^0. (4.38) 
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The second relation in fl4.29l) implies that H^"^^ = H^'^'^ ^ on A, (eiven that 

= there) and hence n{s) is symplectic for Re a: K{s)*{Rea) = Re a. Notice 
that fl4.38l) again shows that is an I-Lagrangian manifold. 

A priori, K,{s){p) is well-defined only for s > small enough depending on p, but 
we shall next derive symbol estimates for ip^^g and that will imply that k(s)(p) 
is indeed well-defined for s small enough independently of p. Assume that^we work 
near a point {xo,r^Q), xq G R", ^ 5"""^, r > and replace (x, ^) by (x, where 

x = Xo + V^x, ^ = r^o+ {r)V^l (4.39) 

Define ^ = A,xo,r,^o by Ai^, = (^i'ix, 0, so that 

i^{x, = ^M^o + v^^, ^^0 + (r)v^O, (4.40) 

d^d^ii =0{l),\/a,(3, for \{x,0\ < Const. (4.41) 
In the new coordinates A, takes the form 



lm{x,0 = ^H^{Re{x,0). (4.42) 



The scaling and the construction of ip^^s commute so if we put 

V'e.s = -^e,s{xo + y/ex,r^o + {r)VeO, 
e 

then f l4.4ip still holds. In particular i/j^^s satisfies the same estimates on A^ as ■j/'e in 
(I4.27P, now with d^, replaced by {dRex, dimx), (c?Re?, <9[m?)- 

The transformation k,{s) can be scaled similarly and the scaling commutes with 
(I4.38P up to a factor (r): If we put 

Kis)ip) = (x(s),e(s)) = ixo + V~ex{s),r^o+{r)V~e^{s)), 
then from fICTD . flCTl) we get 

dsxis) = -^d^ (xis),as)), (4.43) 
dsl{s) = — ^<95;^,,,(x(s),^(s)). 
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We conclude that with (x(0),^(0)) = {y,ri) = (xq + -y/ey, t^o + (f^) y/^V) , 

d',dldl(x,0=Oi{r)-'), (4.44) 
and hence for (x,^) = K{s){y,ri): 

d^d^d^x = C»(e^-^(ry)-'=-l/3|)^ (4.45) 

when k + \a + j3\ > 1. 

Notice that the right hand sides in (14.431) reduce to i{r)~^d^ and —i{r)~^dxip 
when s = 0, where the derivatives are taken in the real directions. The flow in (I4.45P 
is therefore tangent to the one in (I4.32p at s = and we get 

a,XC = 0(e^~^5^(r/)-^-l^l). 
From (I4.45P and the subsequent remark, we have for k{6){x,^) = k{6,x,C,). 

niS, X, = (x, + iSH^, {x, + eU\d{^), (4.46) 

Recalling how k{5) was constructed we conclude that 

k{6, X, = + ^SH^. , (4.47) 

where (x, ^) is real and (x, ^) — {x, ^) = e^5'^{0{-^), 0{-^)). Put as{x,^) = (x, 0- 
The essential part of the discussion above took part near the points of C. In that 
region the discussion is the same in the case when M is compact. 

5 Quantization of weights 

We will follow [13], [23] with one modification; instead of analyticity we will use 
that our weights are "moderate" allowing us to use almost holomorphic extensions. 
Another minor difference is that we shall not use FBI-transforms explicitly, but 
rather rely on certain Fourier integral operators with complex phase. We will assume 
that M = R" for simplicity, but as in the preceding section the essential part of the 
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work will take place near C and here there is no difference between the case M = R" 
and the case when M is compact. 

As a first step towards introducing some Fourier integral operators we shall study 
the function h{y,rj) on T*R", given by 

K{6y{^-dx)-r]-dy = dh. (5.1) 

Recall here that 

= H^^J^{s)p), peAo = T*R^ (5.2) 

where ipe,s is given by fl4.34p . satisfying fl4.33p . Of course, (15.21) remains unchanged 
if we replace ip^^s by an almost holomorphic extension from = A^^^. Now using 
Cartan's formula, we get 

K{s)*{^-dx) = K,{s)*CiH~ ■ dx) 



ds 



iK{sr{H^Jd{i.dx) + d{H^Ji.dx)) 



i^{snH^Ja + d{i.^)) 



Thus, we can take 



m(s)*(d(^ ■ ~ dtlJe,s) 



h = i I (^ ■ ■W7:i'e,s - V'e.s) o n{s)ds. (5.3) 



On As we have = and we recall that ^^{x,$,) = '?/'e(Re (x, ^)). Using (14.341) . 
we get on A^: 



C-Q^ks = ^ ^^e(a;,0 -^^K(Re(a;,e))|(^.)^'50 +^5(«€(Re(a;,e))|(^.):5e), 

(5.4) 

where again we use the notation (■)^, (■)"^ etc to indicate derivatives in the real 
directions. 

We now estimate the terms in the right hand side of (15. 4p . On A^ we have 
Im^ = s{ipey^, so 
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Hence by flOTD . fOH]) . 
Next, look at 

The last term is (5(e(0"^) by dOT]) . The first term is equal to 

(a^K^J^'^ReO = a. ■ g^ii^e)'^ " Re^) - (V^e)'^ ■ a,- 
From flT28|l . fOTD . fOTI) this is d{e{^)-^), so 

(«x|(^.)^'50 = ^(e(ReO-')- (5.6) 

Similarly, 

{a^m"^0 = {a^\{A)^^ReO-^s{a^m'kiAy.) 



-2n 



SO _ 

(«dM?0 = (^(6(Ree)-^). (5.7) 

Returning to fl5.4p . we get 

e-^V^,,, = a(e(ReO"')onA,. (5.8) 

Also recall that ip^^s = V'e so 

e ■ ^^,,s - i^e,s = d{e) on A,. (5.9) 



Combining this with (14. 45 p . we get 



and finally from (15. 3p we obtain 



{^■-^i^,,s-A,s)oK{s) = 0{e), (5.10) 
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Lemma 5.1 

h = d{5e). (5.11) 
Following [13] we shall now quantize k{6) by means of a Fourier integral operator 

Au{x;h)= // e^'^^'^''^'"^a{x,y,a;h)x{x,y,ax)u{y)dyda, (5.12) 

J J (j/,o)eR"xT*A/ 

with a G S^'^"^'^'^'^, and where x is a standard cutoff to a neighborhood of the 
diagonal: |x — a^rl, |y — Q!3;| < 1/C, equal to 1 on a smaller neighborhood of the same 
type. Here we take with P = k{6, a) = K{6){a): 

0(x, y, a) = {x- (3^) ■ (3^ + (a^ -y)-a^ + ^{x - p^,y - a^, a) + h{a), (5.13) 

where on the real domain, 

^l^{x -P,,y- a., a) = 0((ag)((x - + {y - a,f), (5.14) 

in the symbol sense: 

d'.dldZd'^Mx,y,a) = Oi{a^y-\^\i\x\ + |y|)(2-l'=+^l). 



and we take almost holomorphic extensions satisfying the same estimates in the 
complex domain. Further, a G S^'^™''^'^'' means that a = (9(/i~^~'"(a^)t+'=) in 
the symbol sense, dl ydP^d^^a{x,y, a; h) = 0(/i~^~'"(a^)t+'=-l'?l). 
We also assume that 

IrmP^ {a^)iix - (3^f + {y - a^f). (5.15) 

Viewing /i as a function on the graph of k(5), we have dh = f3^ ■ d(3x — «g ■ da^, 
so for 

X = k{6, a)^, y = a^, (5.16) 

we have da4> = 0. Moreover, in a neighborhood of that set, we have with ip = 

i>{x,y,a), 

da(p = (x - K(a)x.) ■ da{K(a)^) + (a^. - y) ■ da^ (5.17) 

dip dip dip 

- — (x - K{a)^, y - a^, a) ■ da{K{a)^) - — ■ da^ + — ■ da. 
ox ay oa 

So is a non-degenerate phase function in the sense of Hormander (a part from the 
homogeneity condition in the fiber variables) with a critical set f l5.16p . the associated 
canonical transformation is k = k{6). 
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Similarly, to k{6) ^ we can associate 

Bv{x; h) = [[ e^?(^'^'^)6(x, y, 7; h)x{x, y, ^)v{y)dyd^, (5.18) 

b G 5'^+"-'?+'', m, /c G R, where with /3 = k{6){^): 

(i>{x, y, 7) = (a; - ix) ■ 75 + (/^x -y)- f3ti + i^{x --fx,y - Px, 7) - ^(7) (5.19) 

and -0 satisfies (15.141) . Again this is a non-degenerate phase and the critical set is 
given by 

x = 7x, y = Px,{P = <S)i^)), (5.20) 

and the associated canonical transformation is k{6)^^. 

For x,y E R", |a; — a^;!, \y — Px\ < 1/C, we have with P = k(q;) = k{6, a): 

lm(j){x,y,a) = (x — Re/5^) Jm/?^ — Im/3a;-Re/5g + Im-?/'(x — /3a;, — a^, /3, a) + Im /i(a). 

Here 

I,n&.Reft = .Re/J,.^ = S(^), 

and 

lmh{a) = 0{Se), 



Imipix - (3^,y- a^,a) > {ai:){— — ^JM_ + _ _ dim (3, 



X\ )l 



so 



Im0(x, y, a) > ^{\x - Re/?,!^ + \y- - C^^^ - C|x - Re/3,;|<5e^ - C5e, 

> ^(|a: - Re/?,|2 + |y - a,n - C5e. (5.21) 
Recalling that e = A/i, A 1, we get 

\ei^i-^y'-)\ < exp (-^(k - Re/3.1' + |y - a.H + C5A). (5.22) 

A similar estimate holds for 0. It follows that A, B are well-defined operators: 
S ^ S with semi-norm estimates that are uniform in powers of h. Also for every 
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s G R, there is an s' e R such that A, B : H' ^ H' with norms bounded by 
some power of h. Moreover, our operators are independent of the choice of almost 
holomorphic extensions of the phase and amphtude modulo operators whose integral 
kernels are (9a (^°°) with all their derivatives and supported in a domain of the form 
\x-y\< 0{l). 

Proposition 5.2 If m,k,m,k = 0, then A,B = C»(exp : L2(R") 

L2(R-) 

The proof will be given later. 

Proposition 5.3 We have BA = Opf^{c), where c = 0{h~'^'^ {a^)^^^) 

Proof 

We have 




where the cutoffs x, X have been incorporated in a, b. Here 

z,7) + = (5.23) 

{x - 7x) ■ 7? + {i^{l)x - z) ■ k(7)^ + (2; - K(a)x) ■ 

+ - y) ■ + '^[x - 7^., z - K(7):r, 7) + ^(2; - y - oi^-, «) + h[pi) - Kl)- 

The contribution to the distribution kernel from a region with |7^ — a^l/ (a^) > 
is 0{h°°), as can be seen by integration by parts with respect to z. More 
precisely the distribution kernel of BA is 

{BA){x,y) = jjj eT^^^+%ax(P^^^)dadzd-i + R{x,y]h), 

where x ^ C^{\Q,oo\) is equal to 1 near 0, and d^dyR = 0{h°°). Here 05.221) and 

the analogous estimate for e^**^ are essential of course. 

Next using fl5.2ip and the similar estimate for Im0, we get a localization in 
|7a: - a^l, leading to 

{BA){x,y) = jjj e^^^+%ax{{l.-a^? )dadzd^ + R{x,y-h), (5.24) 

where d^dyR = 0{h°°). Here x ^ C*^([0, eo[), x = 1 i^sar 0, and eo can be any fixed 
number. 
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In the integral (15.241) . we may assume that |x — \y — ax\, [2; — ^(a)^;! < 1/C, 
where C is as large as we like, since the integral in the complementary region is 
exponentially small. We now want to eliminate integration variables by means of 
the method of stationary phase, and we start by carrying out the z-integration, so 
we first look for the critical point of 

2; z,7) + </)(z,?/,a), (5.25) 

where 0, also denote almost holomorphic extensions. Let Fs{x, 7, y, a) denote the 
corresponding critical value. 

In order to understand this function, we first treat 7 = ^(7), a = K,{a) as 
independent variables (writing K{a) instead of K,{6,a) for short). Let 

G(x,7,7,?/, a,5) = vcz[{x - -f^) ■ -f^ + - z) -^^ + {z - a^) ■ + {a^ - y) ■ 

+V^(x -7x,2;-7x,7) + ip{z -ax,y - a,„a)]. 

Here the critical point z = Zc{x, 7, 7, y, a, a) satisfies 

Zc = ar, + 0{a^ - %, 7^ , {x - 7^)^ {y - a^f) 

in the natural symbol sense. Notice that 

G(7x,7,7; ax,a,7) = 0. 

Moreover, 

dG = 'j^ ■ dx — ■ dy + 'j^ ■ drjx — li • d-^x — (5^ ■ dax — ■ da^) 

^ ^ 7f — Of 
+C((75)(x - 7x, 1/ - ttx, Ix - ttx, ^^^^ )) 

in the same symbol sense, with the convention that the remainder term is expressed 
as linear combination of the "normahzed" forms dx, dy, dax, ctx, d'jx, 7x, {'ls,)~^da^, 



Now, 



and we get 



(x — K{a) 



dFs = -f^-dx-a^-dy + 0{{-f^){x--fx,y-ax,jx-ax,—, — :—)) ■ (5.26) 



a — K,{a.) 
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Moreover, Fs{'jx, 1] 7x, l) = 0. From (14.451) . we know that = ax + 0{5e^^'^ (a^) 
= + 0{5e^^'^) when 3 = K{a), so 

and similarly for "^^^~y^- 
Hence, 

dFs = '^^■dx-a^-dy+0{{'y^)){x-'yx,y-<yx,lx-Oix, ^^ )+d{6){-^x-ax, -^V^)' 

(7?) (75) 

(5.27) 

and integrating this, we get 

Fs{x,'y;y,a) = 75 ■ (x - 7^) - ■ (?/ - a^) (5.28) 
+ (0((75)) + d{6)){x - 7., y - 7. - ^^^(^)'' 

where the loss of e^/^ for each differentiation appears in the variables a, 7 only. When 
5 = 0, we have on the real domain 

ImFs - {a^)iix - 7.)' + (z/ - a^f + (7. - + (^^^(^)')' (^'^S) 

and in view of (15.281) this persists for < 5 ^ 1. 

When applying stationary phase to (15.241) we also have to make a deformation 
of the integration contour in order to pass through the critical point z^. Here we 
recall from [20] and (I5.2ip and its analogue for (p that 



so the error from dz of the almost holomorphic extension, appearing in Stokes' 
formula, is 0{{6e + ImFs)/ (a^))°°). Since e = 0{h), we conclude that 

{BA){x,y) = ||e^^^(™")rf(x,7,y,a;/^M(7.-«x)' + (^^^^ 

(5.30) 

where R has the same properties as R in (15.241) and 

= 0(/i-"^-™-='"+5 (a^)'=+^+^-l^le-^). (5.31) 
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We now compute the Weyl symbol c of BA by means of the formula 

c(x,^;h) = J{BA){x+^,x-^)e-''"^/''dw. (5.32) 

The contribution from R in (15.301) is (9(/i°°(^)^°°) with all its derivatives. The 
contribution from the integral in (I5.30p is 

III e^^^^^^-^^---^^'^^--Mx+^,^^ 

The contribution from a region {\C, — a^ I > ^(a^)} is C»(/i~(0"°° with all its deriva- 
tives and the remaining region can be treated with the method of stationary phase 
by working in the dilated tilde variables given by ■jx = + y/ejx, ■■■ as in the 
addendum below. The proposition follows. □ 

Addendum: Stationary phase with O-symbols. Assume 
e C°°(R"), Im0 > 0, 0(0) = 0, 

|(/)'(a:;)| ^ (j)' = x ■ 0(1), det0"(O) 7^ uniformly in e. 
Let a = 0{1). We shall establish a stationary phase development for 

I{h) = h-^ [ ei'^^''^a{x)dx 



in powers of h = h/e. (Assume /i ^ e ^ 1.) 
Put X = e^x, 

I{h) = ( I ei^^^^a{x)dx, (5.33) 

where (p^e^x) = e0(x), 0(x) = -0(e^x), a(x) = a(e^x). 
Then a = 0{1) in the symbol sense; d"'a = 0{1), and 

\^'{x)\ = \^^'(^ex)\ ^ \x\, ^'{x)=x-Oil) 

and 0"(O) = 0"(O). The contribution to I{h) from \x\ > 1 is 0{h'^) by repeated 
integrations by parts, and the contribution from |x| < 1 can be handled in the usual 
way since = (9(1) here. Thus 

m-p^. c..-g;^e.f-''<»). (5.34) 
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□ 

It follows from the proof that the proposition remains valid if we relax the symbol 
condition in y and only assume 

d'A^ydl^ci = C(/i-^-"^(ag)t+*^-IPle^^(l«l+l^'l)), (5.35) 

i.e. we also allow for a loss of e^/^ for each ^/-derivation. Similarly for B (cf. (I5.18p ) 
we can content ourselves with 

5X,-^74^ = 0(/i-^-'^(«g)5+^-l^le-^(l'?l+lfl)). (5.36) 

Moreover, if 6 an a are elliptic, then c is elliptic. 

We get by standard arguments. 

Proposition 5.4 Let A he an elliptic Fourier integral operator of order (m, k) with 
symbol as in fl5.35p . Then there exists an elliptic Fourier integral operator B of 
order (— m, —k) with symbol as in (15.361) . such that 

BA = 1 + R, (5.37) 

where R is 1-negligible in the sense that its symbol R is 0{{-)'^ {^)^°°) . In particular 
A has the left inverse (1 + R)~^B when e/h ^ 1. 

We notice that when 5 = 0, then A, B are elliptic pseudodifferential operators 
and hence (1 + R)~^B is also a left inverse. (By the Beals lemma we also know 
that (1 + R)~^ is an /i-pseudodifferential operator with symbol 1 + 0{{-^)°°).) For 
general small 6, TZ{As) is closed. Using suitable deformations of A we can produce 
a continuous deformation of closed subspaces in from to TZ{As). All the 
deformed subspaces then have to be equal to and 7l{As) = L^, so (1 + Rs)~^Bs 
is also a left inverse of A^. 

We next turn to Egorov's theorem and start with some preparations. Recall that 
by flOa 

K{a)^ = a^ + 0{5e^ (a?)"^), /«(«)? = a? + 0{5e^. (5.38) 

It follows that 

da{n{(y)f) = dac + 0{6)dax + 0{- — -)dac (5.39) 
da{t^{oi)x) = da^ + 0{- — Ada^ + 0{- — r7:)dai:. 
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Substituting this into (I5.17p . we get 

da(p = {{x - - {y - ar,)) ■ da^ (5.40) 

dip dip 

+ - i^{oi)x,y - Oix, Old ■ 

ox ay 

S 

+ (x - y - a^) ■ {0{S)da^ + 0{- — r)(iag) 

dip, , . , , 

+— (x - y - a^, a) ■ da, 

oa 

which we can write 

da,<P J \{x-K{a),)-{y~a,) + -§^J^''\ad'\ y - a, 

(5.41) 

Here we write a{x,y,a) = 0{m) if d^yd^^dP^a = 0(me-^(lMl+|p|)(a^)-|p|). The 
differential of this vector at a point of the critical set is given by the matrix 

+ <J -{(^d~\€y + ^ + a(^) (5 42) 

1 -1 y (ag) ■ 

If ^ ^ is kernel of the first term, we get = ty and {i^xx+'^yx~^'^xy~^'^yy)'tx = 

0. Here we recognize the Hessian of ^ h- >■ ip{z, z, a) which is invertible because of the 
assumption on the imaginary part. More precisely, the matrix (15.421) has a uniformly 
bounded inverse. From (15.411) . we get 

'^y^ax^)=Msix,y,a) (%^^^"^^) , \x - niaU \y - a^] < 1/0(1), (5.43) 
where 

^ (5 

Msix,y,a) = Moix,y,a) + Oi^), (5.44) 
and d^y^^XMo = Oiia^-W). 

Lemma 5.5 Let ^{x,y,a) = (9((ag)'") in the sense that 
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Similarly, let a = 0((q;^)") be elliptic. Then 

eh'^^(x, y, a)a(x, y, a; h)x{x, y, a)u{y)dyda (5.45) 
eh'^^{K{a)x, oinc, a)a{x, y, a; h)x{x, y, a)u{y)dyda 



eh't'0{\a^)"' ^)a(x, a; h)x{x, y, a)u{y)dyda, 

where x a similar cutoff. 
Proof We have 

The contribution from the remainders to the left hand side in (I5.45P is therefore 

hJJ ^(e^"^) ■ d{e-^a^)^-^)axudyda 
+hJJ ^(e^*) ■ d{e~^ai:)^)axudyda, 

and it suffices to integrate by parts. □ 

Actually, we shall not use the lemma directly, only its proof The next result 
is closely related and could probably be obtained from Lemma 15. 5[ We will give a 
different proof however. 

Lemma 5.6 Under the same assumptions as in the preceding lemma, we have 

eT^'^^axu{y)dyda = AQ, (5.46) 
where Q is an h-pseudodifferential operator with symbol 

Q = $(«:(x,0.,x,x,0 + Oi^iO"^-'). (5.47) 

Here 

Au{x) = eT^'^axu{y)dyda. (5.48) 
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Proof In order to harmonize with Proposition 15. 3[ we may change the assumption 
on a to a = 0{{a^)"^~^2h~~) and assume that a is elhptic in this class. Here O 
indicates a loss of e^/^ for each differentiation in a,y. Let B be of the form (I5.18P 
with b = 0{{a^)~'^^^h ^) where the e^/^ loss is now for each differentiation in 
(x, a). We also assume that h is elliptic. Then from Proposition 15.31 and the remark 
following its proof, we know that BA = Opi^{c), where c = 0{1) is elliptic. Recalling 
that the proof was by stationary phase, we see that if A is the operator given by the 
left hand side of (15.461) . then BA = Opf^(c), where 

c=$(«:(x,o,x,x,oc+a(-(er-'). 

e 

By pseudodifferential calculus, we get BA = BAQ, with Q as in the lemma. More- 
over, B is invertible. □ 



Remark 5.7 Later on, we shall meet the special situation when $ = <!>(«) = 
d{{a^)"'')P{a) =: RP, P{a) = In this case we have 

In fact (anticipating on a part of the proof of Proposition 15. 9p . let 

Bu = JJ eh'^Raxu{y)dyda. 
Then (with Op = Op^ when nothing else is specified) 

5 o Op(P)n = jj ^P{y, hDy) {ei'^Rax)u{y)dyda 

(Piy, -(l>')Rax + d{h^{a^r^^'-^-^))e^'^u{y)dyda. 



As we shall see in the proof of Proposition 15.91 
P{y^ -(i)'y)Raxe^^u{y)dyda = 
jj e^^^axu{y)dyda + jj e^'^d{^{a^)"''+"^'-^)u{y)dyda. 
Applying Lemma [5.61 to B and the various remainder terms, we get (15.491) . 
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Proposition 5.8 Let Ag be a Fourier integral operator quantizing k,{s) as in Propo- 
sition 15.31 Also, assume that a = ag is elliptic and depends smoothly on s in the 
sense that 

where O indicates a loss of e^/'^ for differentiations in a,y. Then 

HDsAs + iAsips = 0, < s < 5o, < 5o < 1, (5.50) 
where Dg = i~^ds and tps is an h-pseudodifferential operator with symbol 

ijs = d{h + h{0-^)+^,^sOK{s). (5.51) 
The term "h" in the remainder can be dropped if ag is independent of s. 
Proof Assume for simplicity that is independent of s. Then, 

hDsAsu{x) =11 eT^'^ads(f)udyda. (5.52) 



Recalhng that, 

(j)s{x,y,a) = {x - k{s, a)x) ■ n{s, a)^ + {ar, -y) ■ a^ + tplx - K,{s, a)x, y - a^, a) + hs{a), 
we get 

ds(f)s{x, y, a) = {x - k{s, a)^) ■ dsK{s, a)^ - dsK{s, a)^ ■ k(s, a)^ (5.53) 

-{dxij){x - k{s, a)x, y - a^, a) ■ dsK{s, a)^ + dshs{a). 

The restriction to the critical set (15.161) is 

(9s0s)(ft;(s,a)x,aa;,a) = {dshs){a) - dsK{s,a):, ■ K{s,a)^ (5.54) 

d - 

= {dshs) {a) - ^i-g^i'^,s ■ ° '^(■5) 

= -i^lJe,s ° K,{s){a) 
= 0{e), 

where we used ( 15. 3p . More generally, from (15.531) we get 

ds(f)s = -i^e,s o K{s){a) + d{e^) ■ {x - K{s,a)x,y - a^;)- (5.55) 
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As in the proof of Lemma (15.51) we can make integrations by parts and see that the 
contribution from the remainder to (15.521) becomes 

j j ei'l'd{h-"'-^{a^)''+^h{ai:)-^)u{y)dyda. 

Combining this with (I5.55P and Lemma 15.61 we get the proposition. □ 

We can now prove Proposition 15. 2[ 
Proof We only consider A and we may assume that A = A{6) where A{s) is a 
smooth family as in the preceding proposition. The result for B will be the same 
since B is like the adjoint of A. From (15.501) we get 

hd,A: = rsA:. 

If M G iS, then A*u G S and we get 

hds\\Aluf = {ijlAlu\Alu) + (v4*m|^*A», 

so 

hds\\A*uf < 2||V^:||||v4>||^ 

But 11^:11 = C(e), so 

\\A:u\\ < e^(^)l^l/'^||w|| = e^«^l^l||w|| (e = Ah). 

Then \\Ag\\ < e'^(^)"^'^ and the Proposition follows since As and A} have the same 
norm. □ 

Let P = in the symbol sense: d^d^P = We study PAs - 

AgP where Ag is as in Proposition 15. 8[ Sfa^ bounded in S™'~^^'^~^^ (so no -^/e-loss 
in the symbol for simplicity). Working with the Weyl quantization we know that 

P{x,hD^;h){ei^'^^'y'''^a{x,y,a;h)) = + d{W^ {a^)-^) 

'P{y,hDy;h){e^''''^'''y'"'>a{x,y,a;h)) = e^<^c + 

where 

b = P{x,<P'^{x,y,ay,h)a{x,y,a-,h) + d{h'''^~'''{a^)'+^+^'h\a^)-^), 
c ^ P(y,-0;(x,i/,«);/i)a(x,i/,«;/i) + a(/i-^-™(«5)'=+t+A^/i2(c,^)-2)^ 

and the O refers to -^e-loss only with respect to differentiations in a. (Recall the 
general fact that the Weyl symbol of e~"^' o p'^[x, hDr) o e*"^/^ is equal to p{x, ^ + 



36 



We conclude that PA, - A,P = A, + A, o Op {0{h'^{ai:)-'^)), where 

A,u= 11 ei'^{P{x, (j)'^{x, y, a);h)- P{y, ; h))au{y)dyda. (5.56) 



On the critical set fl5.16p . we have 

P{x, 0^; h) - Piy, h) = P{n{s){a)) - P{a), 

and more generally, 

P(a;,0;)-P(y,-0;) = P{K{s){a)) - P{a) 

+ {x - k{s, a),) ■ {0{{a^f') + d{se^ {a^)''-')) 

Using flOal) . flOil) . we get 

P(x,0;)-P(y,-0;) = PKs,a))-P(«) (5.57) 

The contribution from P(/t(s)(a)) — P(«) in fl5.56p is of the form 

A, o Op {Pi^{s){a)) - P{a) + d{s^{a^Y'-')) 

in general, and the remainder estimate improves to 0{sh{a^)^~'^) if P{K{s,a)) — 
P(a) = C(se(a)^~^). By integration by parts in a, we see that the contribution 
from the two remainder terms in fl5.57p is 

hjj e**5„^((0((a5)^) + a(se^(«^)^-i)+9„,(^(idem))a)w(y)d2/da 
= h [[ ei'^iOiiai:)^'-') + d{s{a^y'-^))au{y)dyda. 



In conclusion. 

Proposition 5.9 We have 

PAg = AgPg, Pg = P -\- Qs, 

where 

Qs = P{K{s){a)) - P{a) + a(s4(a«)"'-') + d{h{a^)''-') 

£2 
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in the general case, when P(/€(s)(a)) - P{a) = 0{se2 (a^)^^-^) . 

In the special case when P{K{s,a)) — P{a) = 0{se{a^)^^~^) the first remainder 
term improves to 0{sh{a^)^^~'^) so in that case, 

Qs{a) = P{k{s, a)) - P{a) + d{h{a^y'-'). 

In order to treat certain conjugations, we need a more precise description of Qs 
in the general case in the last proposition. The proof above shows that 

AsQsU = ff ei'^{P{K{s, a)) - P{a))au{y)dyda + A,Op (5.58) 



and we need to take a closer look at the oscillatory integral. 
By Taylor's formula and (14. 45 p . 

P{K{a)) - P{a) = OisyTe) ■ (P;(a), -^fLi^)) + d{sh{a^)''-'). 
When passing to Weyl composition of symbols, we notice that if r = 0{s^/e), then 

Lemma 15.61 and Remark 15.71 then show that on the symbol level 

Qs = dm)''"') + disMO''-') + 0{sVimP'^. ^P^). (5.59) 

Here the first term to the right is too large; we would like to have h{^)'^^~^, so we 
take a closer look at Qs using Proposition 15. 8[ where we now add the assumption 

m = k = 0, Aq = 1, is independent of s. (5.60) 

Then as noticed, (15.511) improves to 

i^s = i^sOK{s)+di^) = die). (5.61) 

We get on the symbol level, writing [A, B] = Ajj^B — Bjj^A: 

hdsQs = [Ps,i's] = [P,^s] + [Qs,^s] 

= [p, iJs] + dm^-h^) + d{sh{ir^"h^) 

+(5(.v^e-|y)#(P^, ^P^) + a(.v^)#a(e(o"^-^)^) 



[P,^s] + O{h'{0^''-') + O{s'eh{0 
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On the other hand, 

[PM = -{P,^s} + d{{0''e 



thiii, o k{s), P} + hO{^^mP'^, ^) + 0(^(0"'-') 



so 



1/2;^ pi 

hdsQs = ih{^^sOK{s),P}+0{-^^+^mP'^,jl^)+0{h^i)''-^ 
Since k{s) is obtained from integrating zif^ , we obtain, using also that Qq = 0, 

= PoK{5)-p+oC-^+'-^mp^, f^)+omO''-'+sMO''-') (5.62) 

under the assumption fl5.60p . (Recall that here M is the order of P.) 

6 The conjugated pseudodifFerential operator 

Let P be the operator introduced in the beginning of Section [2] satisfying the as- 
sumptions fl2.1l) - fl2.7p . To start with, we assume Hypothesis 14.11 and define C,'?/'^ 
as in (14.51) . (14.61) . Later, we shall also use the dynamical Hypothesis 14.31 (implying 

MD)- 

The function satisfies V. = 0{e), ^ ■ d^tp, = C(e(0"^) in the sense of f OTD . 
(I4.28p . but we also know that these bounds improve considerably away from C. 

Let A = As he an elliptic Fourier integral operator of order 0,0 quantizing k{6) 
as in Section [5l and let Bs be a corresponding operator quantizing k{6)~^. Then we 
can apply Proposition 15.31 to see that BA = Opf^{c), where c = C(l) is elliptic. Now 
restrict the attention to the intermediate region, where ip^ = etplx), iplx) G C°°. 
Using that e = Ah with A large but fixed, we find here that 

where c = 0{1) in the standard symbol sense, d^d^c = O((0""^') and c is elliptic 
in this class. Let c = 0{1) be an elliptic symbol which is equal to 1 in the interior 
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region and which changes from 1 to c when going outwards in the intermediate 
region. We now replace Ahj Ao Opf^(c)~^. (Here we assume that Aq = 1, to avoid 
a topological difficulty.) Then we have achieved that Au = e^'^^/^u for u supported 
in the exterior region including a shell where ip^ = eip , ip = il){x) G as above and 
the region further out where ipe = ipe{,x) = 0{e). 

Next, we consider defined as in Proposition ED by PAs = AsP^. Recall that 
P = P2 + tPi + Po and write P^As = A^P^ , P^ = P + Q^, P/ = Pj + Qj. 

For Pj, j = 0, 2, we apply fl5.62p and get on the symbol level over any bounded 
set in x-space 

= diSV~emd^p„^) + diihS + Sh){0'-'). (6.1) 

For pi we need the more precise information about ip^ given in (14.101) . (14.120 . 
(I4.13p . satisfied also by ipe, as well as the fact that ipe{p) is independent of e in 
\p~ Pj\ ^ (actually in \p — Pj\ < \^/C for some C > 0, but we can always dilate 
in e). From these estimates it follows that 

%i) ® {d.A, iOd^A) = d{e) (6.2) 

over a neighborhood of the set where ip^ also depends on ^. It also follows that 

Pi(«(5)(x,0)-Pi(x,e) = 5(5e). (6.3) 
From Proposition 15.91 we deduce that 

Qi=p,{Ki5)ix,0)-Pi{x,0 + Oih). (6.4) 
By construction of k{S), especially (14.471) . we see that 

P,{k{6){x,0)-Pi{x,0 = iSH^M^.O + 0{e5''). (6.5) 

In fact, this is quite obvious in the intermediate and far out regions, and near a 
point Pj G C, we choose canonical coordinates so that pj = (0, 0) and put p = \fep. 
Then if '?/'e(p) = e^e(p)5 '^^ get = H^^, where the tilde on the H indicates 
that we take the Hamilton field with respect to the p-variables. The manifold As : 
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Im p = 6H^^(Re p) now becomes Im p = 6H^^{Re p) while pi = 0{e]?) and since tp^ 
satisfies fHTTUl) . 

Then k{5){p) becomes k{S){p) = p + iSH^^{p) + 0{5^{p)-^), so Pi{k{5)p)) = p^{p) + 
iSH^^ipi) + C(e52), leading to (lO) . 
It follows that 

Qi = -t6Hp^^, + d{6h + h) = d{5t + h). (6.6) 

Away from any neighborhood C we have the improved estimates 

Q{ = -i6Hp^ij, + 0{6h + h) 

in the usual symbol sense, as long as we stay away from the outer region where 
ipe = ipe{.x) only satisfies -i/^e = C(e), and where we can apply the analysis of Section 

m, so 

Q{ = -i6Hp,ij, = -iSuiij,) (6.7) 

there. 

7 Estimates for the conjugated pseudo differential 
operator and localization of the spectrum 

Let be the conjugated operator of Section [61 We shall study lower bounds for 

Re {P^u\u) = Re {{P^ + P2)u\u) - Im {Piu\u), u e S{M). (7.1) 
Using (16.11) . we get 

Re{P^u\u) > (P2w|w)-C5v^||w||||Op((a^P2,^))w||-C(M + (52e)||wf 
> iP,u\u) - ^||0p ((%2, ^))ur - C2{h6 + 6h)\\ur, 
where Ci can be chosen arbitrarily large and C2 depends on Ci. Here 

2P2-^(5,p„^) >0, 
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if Ci is large enough, so 

{P2u\u) -^j\op mp,,^))ur = (op(P2-^(%2,^)>w+0(/i^)hi 

> ^iP2u\u)-Oih')\\ur, 

where the estimate follows from the Fefferman-Phong inequality in the semi-classical 
setting (see [Rj). It follows that 

Re{Piu\u) > ^{P2u\u)-C{h'^ + h5 + 6h)\\uf. (7.2) 

For Pq we have the same conclusion, 

Re iP^u\u) > ^{Pou\u) - C{h^ + h6 + ^^e)!!^!^, (7.3) 

since the same analysis applies for x in a bounded region and further out, we just 
have P^ = Po>0. 

Combining this with (16.61) . we get 

Re {P^u\u) > (5(0p {Hp,ij,)u\u) + ^((Pq + P2)u\u) - C{h + h6 + Sh)\\uf. (7.4) 

Here we recall that Hp^ip^ = {pe)To ~ Pe by (14.71) . that (p'e)ro satisfies (14.251) . (I4.26P 
and that Pe < P < Po + P2 i>Y (I4.16p . Write, 

SHp^i:, + ^(po + P2) = HPe)To - Spe + ^{po + P2) (7.5) 

= ^(Pe)To + Kp-Pe) + 5(P0 + P2-p) + {^-5){P0+ P2) ■ 

Here we want a lower bound for Op {p — Pe)- This is quite straight forward away 
from C, so we concentrate on a neighborhood of a point pj G C. Assume pj = for 
simplicity. Then near we have by (14.151) 

|p|2 _ _ _ |p|2 
Pe = 9i — )P, P-Pe = {I- g{ — ))P = (1 - ^e{p)Yp, 

where 



42 



and we may assume that g has been chosen so that 1 — g = {1 — i)"^ with i smooth. 
Here 

(i-4)#P#(i-4) = 

= {l-i,{p)fp + 0{h). 

From this we conclude that 

(Op {p- p,)u\u) > ~Ch\\uf. (7.6) 

Similarly, 

iOpipo+P2~p}u\u) > -Ch^Wuf, (7.7) 

by the Fefferman-Phong inequality or by a direct argument. 

Let < = 0{e) be equal to e in C + -8(0, y/e) and have its support in C + 

B(0, v/2i). Let K, = Op (K). By fOSD . flCT]) we have {p,)to ^ min (e, dist (-, CY), 
so _ 

e^K + {p,)To = Oie). (7.8) 

Hence, 

{{K, + OpmT,))u\u)^e\\uf. (7.9) 
Combining fl7.4l) - fl7.9p . we get 
Proposition 7.1 VFe have 

Re ((P-^ + ir,)M|M) > ^llwf + (^ - 5)((P2 + Po)u\u) - C/i||Mf , M G 5(M). (7.10) 

Here we recall that 6 > should be small enough, e = Ah with A arbitrarily large 
and fixed, and C in (17.101) is independent of S, A while h is small enough depending 
on these two parameters. From fl7.10p we get the a priori estimate 

-Ch- Rez)\\u\\ < \\{P^ + K,- z)u\\, u G 5(M), (7.11) 

when Rez < 5Ah/C - Ch. 

From Section [3] we know that P has no spectrum in the open left half-plane. We 
shall next prove 

Proposition 7.2 For every constant B > there is a constant D > such that P 
has no spectrum in 

{z eC;Rez < Bh, \lmz\> Dh} (7.12) 

when h > is small enough. Moreover \\{P — z)^'^\\ = Osih"^) for z in the set 
(EI2]). 
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Proof Choose S > small, then A large enough, so that ^^ — Ch — Kez > e/Const. 
when He z < Bh. Then (17.111) gives 

T^lkll < ||(^^ + K,- z)u\\, u G S{M). (7.13) 

Take z in the set (17.121) . When Rez < 0, we already know that z ^ o"(-P), so we 
may assume that Q <'Re z < Bh. 

Now recall that the symbol of is 0{e) and supported in C + 5(0, \/2t). On 
that set we have = (9(e) = 0{Ah) and hence 

Ip'' — z\> — — — , when \lmz\ > Dh, 

and D is large enough, assuming still that < Re 2; < Bh. 
It follows that we can find E = 0{t/{t + \z\)) such that 

K, = Eo{P^ ~ z) + F, F = d{h), 

where E, F also denote the corresponding /i-pseudodifferential operators. In partic- 
ular, 

\\K,u\\ < - z)u\\ + Oih)\\u\\, (7.14) 

where the Os are uniform in e. Combining this with (17.131) with e ^ h, we get 

e||?i|| < 0(1)||(P'^ - z)ul u e S{M), (7.15) 

for < Re 2 < Bh, < Imz < Dh. Now recall that = A^'^PAs, where 
As, AJ^ -.S ^ S, E^ and have norm Oa(1). Then (17A5|l gives 

h\\u\\ < Ob{IMP - z)ul u e S{M), (7.16) 

for z in the set (I7.12p . From Section [3] we then know that z ^ cr(P) and that 
\\{P~z)-'\\<OB{l)/h. □ 



8 Asymptotics of eigenvalues 

Let pj G C and let Fp^ be the matrix of the linearization of Hp at pj (the so called 
fundamental matrix of p at the doubly characteristic point pj). Thanks to the fact 
that the quadratic approximation of = p o k{5) at pj is elliptic on Tp.{As) and 
takes its values in a closed angle contained in the union of {0} and the open right 
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half plane, we know from [23] that the eigenvalues of Fp. are of the form ±Aj^fc, 
1 < k < when repeated with their multiplicity, with Im > 0. Let 



tr {p,pj) = -^Aj-fc. 



In our case the subprincipal symbol of P at pj is zero and will not enter into the 
description of the eigenvalues. 
Put 

lix, = -Pi^^ ^0 =P2+Pl- PO- 

Let Fq, Fp be the fundamental matrices of q,p at one of the critical points pj G C. 
Since 

Id d 

we see that Fg and ^Fp have the same eigenvalues; ±jXk, k = l,...,n (j being 
fixed) where Re(iAjt) > 0. Now q is real-valued and we can apply the stable 
manifold theorem as in [10] (and at many other places) to see that the if^-fiow has 
a stable outgoing manifold A+ passing through pj such that TpA^ is spanned by 
the generalized eigenvectors corresponding to +jAa:, A; = 1, n. We also know that 
A+ is a Lagrangian manifold and that q vanishes on A+. 

Lemma 8.1 Assume for simplicity that pj = (0, 0). Then Tp.A^ is transversal both 
to {x = 0} and to = 0}. 

Proof Since we consider the linearized situation we may assume right away that 
g is a quadratic form (the second order Taylor polynomial at pj), so that pj are 
quadratic forms as well. The dynamical condition fl4.2ip implies that 

{po +P2)to > as a quadratic form. (8.2) 

Let L = A+ n {x = 0}. Since po = 0, pi = on {x = 0}, we know that p2 = 
on L and hence Hg = Hp^ on L. Now Hg is tangent to A+ and Hp^ is tangent to 
{x = 0}, so Hg = Hp^ is tangent to L. Thus L is an Hg- and Hp^ invariant subspace 
on which P2 + Pq = 0, and since {p2 + Po)to > away from 0, we necessarily have 
L = 0. 

The proof of the fact that A_|_ fl {.^ = 0} = is the same after permuting the 
roles of Po and p2. □ 

It follows from the lemma that 

A+: ^ = (P'4x), XG neigh (0), (8.3) 
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where 0+ G (neigh (0); M), 0+(O) = 0, 0'+(O) = 0, det0'|(O) ^ 0. 

Let A_ be the stable incoming ifg-invariant manifold such that Tp.A^ is spanned 
by the generalized eigenvectors of Fg corresponding to —jX^, 1 < k < n. The lemma 
is valid also for A_ and fl8.3p has an obvious analogue for A_ where we let 0_ denote 
the corresponding generating function. 

Proposition 8.2 We have 0'|(O) > 0, 0" (0) < 0. 

Proof Again we can consider the linearized quadratic case. If we make a smooth 
deformation of q, then A_|_, A_, 0+(O), 0" (0) vary smoothly with the deformation 
parameter and det0'4(O) ^ 0, provided of course that we maintain the condition 
(18.21) . Consider the deformation from g = go to to ^"^ — x'^ = qi. 

qt{x, = (1 - t)q{x, + tie - = V\ + V\ - pI 

with 

P2(0 = (1 - t)P2{0 + P'oix) = (1 - t)poix) + tx\ 

p\ and Pq are positive definite for t > 0, so (18. 2p is maintained. For t = 1 we have 
0j_(x) = ±a;^/2 so ±0'4(O) is positive definite. Since the signatures of 0'-t(O) are 
independent of t, we get the lemma. □ 

From Section [6] we recall that the conjugated operator Ps = Ps,e has the symbol 
p + SH,,ij, + d{6^emd^p„ ^) + d{h + 

j=0,2 

and that we have the a priori estimate (I7.10p expressing that the real part of 
is > % - Ch outside C + 5(0, ^e). In the set C + S(0, v^e) the symbol P^ is 
independent of e modulo (9(e(^)°°) and is of the form ~ p'' + hrx + /i^r2 + 
where ^ dist(-,C)^, Re p^ - dist (^C)^. In the following we shall assume for 
a while (in order to simphfy the notations) that C is reduced to a single point 
Po = (0,0). 

We fix i? > and restrict the spectral parameter z to the disc -D(0, Bh). We shall 
take e = Ah with B sufficiently large. Assume that M = R" for simplicity. 

Following basically [15] we recall the construction of a well-posed Grushin problem, 
first for Ps — z and then for Ps — z. Let A = (x^ + {hD^YY^'^ so that Ps is equipped 
with the natural domain 'D{Ps) = {u e L^; k^u G L2(R")}. In Section 11 of[15] the 
authors constructed operators 

R_ : c^" R+:L^-^ (8.4) 
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of the form 

No 

R.u_ = 5^w„(j>J(x), RMj) = Nf^i^)), (8.5) 

with the following properties: 

eJ(x) = /^-?e,(-^), /;(x) = /.-t/,(-^), (8.6) 

e,{x) = j9,(x)e^*°(^), f,ix) = g,(x)e^*«(^\ (8.7) 

where pj, qj are polynomials and $0? are quadratic forms with Im $0; > 0, 

$0 ~ If (5 > is small and fixed, A sufficiently large, the problem 

{Ps - z)u + Rtu. = V, R^u = v+, (8.8) 

for V G L^(R"), t>_|_ G has a unique solution u G V{Ps), G C^", where 
i?^ = -R+A^, i?l = y4^^i?_. Moreover, for the solution, we have the a priori estimate 

||A^m|| + \uA < cdl^ll (8.9) 

Notice here that {h-^''^A)^ R. = 0(1) : C^" L^, for every G R and similarly 
for R^. From this, it follows that Rt = R^- depends weakly on e in the sense 
that if e = Ah, A> A, then 

and similarly for R^. 

We shall derive an a priori estimate for the problem 

{Ps - z)u + Rtu^ = V, R^u = v+, (8.10) 

when u E S. Let x ^ ^'^(5(0, 2)) be equal to one on _B(0, 1), and put X^i^^O = 
x(e~^/^(a;, ^)). We use the same notation for the corresponding /i-quantization. We 
may assume that Pg — Ps = C(e(/i/e)°°) on supp (Xy/l)- From the first equation in 
fl830|) . we get 

{Ps - z)x^eU + Rtu^ = X^eV + X^e{Ps " A)m + [^5, X^^U + (1 - X^^R-U-, 

R+X^eU = f + - i?+(l - X^)u. 

(8.11) 
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Here 

11(1 - Xv^)R'-U-\\ = 11(1 - = (8.12) 

h\R'^{l - x^)u\ = /i|i?^(/i-U)^(/^-U)-^(l - xv^)u\ = 0{{^)^h\\u\\). (8.13) 
Thus, applying the a priori estimate 08.91) to fl8.1ip . we get 

\\A\^M\ + < Ci\\xv-M\ + Oih^)\\u\\ + Oi^)\u^\ + h\v+\). (8.14) 

We next look for an a priori estimate for (1 — x^)m. Apply 1 — to fl8.10p : 
{Ps -z){l- x^)u = (1 - Xv^)v + [Ps, Xvi]u - (1 - Xvi)R~u^- (8.15) 
As before, 

[Ps, xv^] = [Ps, XV^] + [Ps - Ps, XV^] = 0{h^) : ^ L\ 
and using also (I8A2|) . we get from (|212|) . (IHTTSD that 

^11(1 - < 11(1 - + \\m - Xv/i)«ll + ^^(^Mhll + Oi^)\u^l (8.16) 

where we can take K^{x,^) = ^x C^^^ )- Since 

we get 

^11(1 - XV^M < 11(1 - XV-M\ + (^(4)lhll + 0{-J\u^\. (8.17) 

Let = 0{mm{h + dist (x, ,^ ; 0, 0)^, e)) be an elliptic symbol in the class defined 
by the right hand side, so that 



\\^\v-M\+4i^-xv^M-\\LM 
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Then summing (I8.14p . (I8.17p . we can absorb the various remainder terms to the 
right, and obtain 

^(||L,Mi| + |n_|) < \\v\\ + h\v+\, (8.18) 

when f lSTTOl) holds. 

Now recall that Ps, Ri have been defined from P, i?_ by conjugation 
with As, and use also that 

^h\\u\\ < \\L,u\\ < CAh\\ul 
to see that if u^v ^ S{M) and 

SiP-z)u + R_u. = v ,^^(0^^;,)^ (8_19) 

V4- 



\r+u 



then 

h\\u\\ + < C{\\v\\ + /i|t;+|). (8.20) 
From the discussion after the proof of Corollary 13.21 we conclude that 

V(z) = {jj^^ ^' ) : T^{P) X C"^" ^ X (8.21) 

is bijective with a bounded inverse 

8 = ^^^^l^ ^^^^l^ y. L^x C^o D(P) X C^o, (8.22) 

for z e /^(O, P/i), and (Km shows that 

||i?(z)|| = 0(i), ||i?_(^)ll = (^(1), \\E4^)\\ = 0{1), \\E_4z)\\ = 0{h). (8.23) 

In Section 11 of [15] the authors studied the action of V{z) on spaces of functions 
of the form {a{x; /i)e**°*^^^/'^, m_), where a is a symbol, and deduced that E_+(z; h) 
has an asymptotic expansion in half powers of h with a certain additional structure. 
From that was obtained the asymptotic expansion of the zeros of det i.e. of the 
eigenvalues of P in D{0,Bh). That discussion goes through without any changes 
in the present situation, so we get the asymptotics for the eigenvalues in any disc 
D{0,Bh), when h ^ 0. 
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Theorem 8.3 We make the assumptions fl23]) -dO) . 1^^, (g^H), i K2^ . and re- 
call the definition of C in f l4.5p . Let B > 0. Then there exists Hq > such that for 
< h < Hq, the spectrum of P in D{0,Bh) is discrete and the eigenvalues are of 
the form 

Xj,kih) ~ hifij,k + h^'^^'-'HAi + h^'^^'-'NA^ + •••), (8.24) 
where the fij^k are all the numbers in D{0, B) of the form 

1 " 1 ~ 

A^i.fc = ~ 5Z ^j'f^'^^j'^ + Pj^' ^^^^ ^i''''^ ^ ^' (^-2^) 

l=\ 

for some j G {1, A^}, = #C. (Possibly after changing B, we may assume that 
7^ B, ^j,k.) Recall here that ±Xj are the eigenvalues of Fp. This description 
also takes into account the multiplicities in the natural way. If the coefficients Vj^k^i 
in fl8.25p are unique, then Nj^k = 1 o,nd we have only integer powers of h in the 
asymptotic expansion (18.241) . 

Theorem 8.4 We make the same assumptions as in Theorem 18. 3[ For every B, 
C > there is a constant D > such that 

\\{z-P)-'^\\ < ^, for z e D{0,Bh) with dist (z, a (P)) > ^. (8.26) 

The last result follows from the formula 

(z - P)-i = ~E{z) + E4z)E_+iz)-'E4z), 

(K23\f and the fact that ||^Z+(-2)|| = 0{h^^) when (18:261) holds. 
Still with j = jo fixed, let 

1 " 1~ 
/i = - 5^ ^^^^ + (P' Jo), e N (8.27) 

be a value as in (18.25!) and assume that fi is simple in the sense that (i^i, i^n) £ 
is uniquely determined by /x. In particular, every for which i/£ 7^ is a simple 
eigenvalue of Fp. Then as in [10] (see also Chapter 3 in [5j) we can construct 

X{h) h{fi + hfii + h'^fi2 + ■:) (8.28) 

with uniquely determined coefficients fii, /i2, ... and 

a{x; h) ~ ao(a;) + hai{x) + ... in C°°(neigh (xjo)), (8.29) 



50 



where aj{x) = 0{\x — XjJ^"^ ^"'''"''), = ^ ^'^ and ao has a non- vanishing Taylor 
polynomial of order m, such that 

(P - A(/i))(a(x; /i)e-^+(")/'^) = C(/i°°)e-*+(")/'^ (8.30) 

in a neighborhood of Xj,y Actually any neighborhood Q CC R" will do, provided 
that 

1) is well-defined in a neighborhood of Q. 

2) H,^^^^OonQ\{x,,}. 

3) fl is star-shaped with respect to the point xj^ and the integral curves of the vector 
field Z/+ := (vr^)^ {Hg, ), where 7i^{{x,^)) = x. 

We also know that X{h) is equal mod 0{h°°) to the corresponding value in (18.241) . 
As in we notice that if 7 C Z)(0,i?) is a closed /i-independent contour 
avoiding all the values /x^ fc in (18.251) . and 

T^h'y = ^ [ {z- P)-'dz (8.31) 

the corresponding spectral projection, then, using also Theorem 18. 4[ 

||7r,^(xae-^+/'^) - x^e-'^+Z'^IU^ = 0(/i°^) (8.32) 

if X £ C'^(f2) is equal to one near xj^. It follows that xae~'^+/'^ is a linear combination 
of generalized eigenf unctions of P with eigenvalues inside hj up to an error 0{h°°) 
in L^-norm. 

9 Exponentially weighted estimates 

In this section we keep the general assumptions on P and assume for simplicity that 
C is reduced to a single point: 

C = {(0,0)}. (9.1) 
Hip e (neigh (0, M); R), we have 

e^/^ o P o e-^/^ = P^, (9.2) 

with the symbol (cf fl27[3D ) 

P'4>i^, = P2ix, - li^i ^'(x)) + iid^(l)ix, ij'ix)) ■ ^, (9.3) 
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where we recall that 

q{x, i) = p2ix, + Piix, - Poix). (9.4) 

Notice that ^ t— q{x,C,) is a convex function for every x. 

Let (j) = (t)+{x) G C°° (neigh (0; R)) be the function introduced in Section [8] so 
that A+ = is the stable outgoing manifold through (0, 0) for the iJ^-fiow. Recall 
that by Proposition 18.21 

0"(O) > 0. (9.5) 

We have the eikonal equation 

g(x,0'(x))=O, (9.6) 

so P4,{x,C,) = P2{x,^) + i{d^q){x, (p'{x)) ■ ^. The vector field {d^q){x, 4>'{x)) ■ dx is the 
x-space projection of Hq^^ , so its linearization at x = has all its eigenvalues with 

real part > 0. Consequently (as we shall see in more detail in the proof of Lemma 
OroU below), there exists G G (neigh (0, M); R) such that 

(%)(x,0'(x)) ■ d^G - x^, G{x) - x^. (9.7) 

Let Vtcir) = {x G neigh (0); G{x) < r} for < r < 1. Outside the set fiG(Coe), 
we put 

ij = <i)-eg{G), (9.8) 

for a suitable function g and for < e ^ 1. (Eventually e will be proportional to 
h.) Using (19. 6p . we get 

g(x, ^\x)) = -tg\G){d^q){x, 0'(x)) ■ 9,G(x) + t^g\Gfp2{x, G\x)). (9.9) 

Choose g{G) = InG for G > G^e, so that g'{G) = 1/G. Then 

,(.,V7(.), = _.J_^(a,,,(.,^'(.), .4G + .^?!H%|M. (9.10, 

Here 

i(%)(x,</)')-5.G-^l, P^^^ = Oil). 



We conclude that 



g(x,?(x)) < xG neigh (0,M)\fiG(Coe), (9.11) 
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if Co > is large enough. 

Outside a small fixed neighborhood of we want to flatten out the weight. Let 
fs{t) = ^/(|) be the function introduced in Section [2l For some small and fixed 
Sq > 0, we put 

i^ = fsM = fsM-^9iG)) (9.12) 

which is also well-defined as the constant 36o/2 for large x. From (19. lip , the fact 
that q{x, 0) < and the convexity of q, we get 

qix,^P'ix))<~fs,W, (9.13) 

where we keep in mind that < /^^ < 1. 

We extend the definition of to a full neighborhood of x = 0, by putting 

g{G) = In(Coe) + -^{G - Coe), for < G < Coe. (9.14) 

Then in ndCoe), we have ^' = = (f)' - G'/Co, so g(x,V^') = 0(e) by ([93]). 

In the exterior region where ip' = fsg{^)4'' is small, the conjugated operator 
P = is close to the unperturbed operator P and we can apply the method of 
Section H Write 

p = P2 + ipi+Po (9.15) 
for the symbol of P, so that by (19.31) . we have 

Pi = {d^(l){x,'ip'{x)) po = -q{x,ij'), (9.16) 

to be compared with pi = {d^q{x,0)) ■ ^, po = —q(x,0). Let iy{x,dx) = J2^jix)dxj, 
where pi = '^Cj{x)^j and define ip^ as in (12.151) . now for P instead of P and with 
an additional cut-off: 

= (1-X(x)) J ki^)f,opooexpitu)dt. (9.17) 

Here x ^ C'^(M) is equal to 1 near x = and has its support in a small neighbor- 
hood of that point. 

Then for Pj^^ = e^^^/'' o P o e"^^^/'* we get (cf (E^Zl)), 

Ps^^ > (1 - S)po + 6{feO Po)to - 0{6^e) outside supp x, (9.18) 
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where the time average in the second term to the right is taken along the trajectories 
of v. In the region where x = 1, "we get = P of course, and in the intermediate 
region, supp (Vx), we have 

PoM.=Po + 0{5e). (9.19) 
If 5i > is small enough, we know that 

/;(^(x))<5l^(/eOpo)To(x)>^, 

for some constant Ci > 0. In fact, the u and u trajectories through a given point 
with fsg{i^{x)) < 6 stay close for some fixed time > , so the conditions (14.241) . 
(I4.23P imply that the //-trajectory will encounter points with Pq > 1/Const during 
a non-trivial interval of time. 
Then 

• In the region where fs^i^ix)) < Si, we have p^g^^ > (1 - S)po + ^ (and we 
recall that po > 0). 

• In the region where /^^ ("^(a;)) > 6i, and G{x) > Cqc, we have Pg^^^ = Po + 
0{6e) >^ + 0{6e), by I^JS^. 

• In ndCoe), we have = Po = 0{e). 

Choosing first 6i > small enough, then 6 > small enough, we conclude that 

J > ^outside naiCoe), 

Now Pg^^ = e^^l^ oPo e-^-l^ = P^^, where 

^, = ^ + (9.21) 

and where we recall that ip also depends on e. Combining Lemma 12.11 for P^^ with 
dnSni), we get 

/ {-^-Rez)\u\'^dx+ / {P2u)udx< (9.22) 

Re{{P^^- z)u\u) + I {0{t) +Rez)\u\^dx. 
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If Re 2; = 0{h), we choose e = Ah with A so large that — Rez > h and deduce 
that 

h\\uf + iP2u\u) < Re ((P^, - z)u\u) + O{h)\\u\\l^^coe). (9.23) 
leading to the a priori estimate 

h\\u\\ < ||(P^^ - z)u\\ + O{h)\\u\\na^coe). (9.24) 
Re-injecting this estimate in (I9.23p . we get 

h^uW" + hWBhDuf < (1 + a)\\{P^^ - z)nf + O^ih^)\\u\\l^^coe), (9-25) 

for every fixed a > 0. 

Here (P^^ — z)u = e'^^^'^{P — z)v, u = e'^^^^v, so from fl9.25p we get the a priori 
estimate for the original operator 

/i||e^^/^|| +/i5||pUD(e^^/''i;)|| < Oil^e'^^^^P - z)v\\ + Oih)\\e^^^''v\\na(Coe), 

(9.26) 

uniformly, for |Rez| < Ch provided that e = Ah for A large enough depending on 
C. 

Now let X{h) = \i^k{h) be an eigenvalue of P as in fl8.24p . fl8.28p and assume 
that n is given by (18.270 and is simple, as explained after that equation. Then X{h) 
is a simple eigenvalue of P and is the only eigenvalue in some disc D{X{h), h/Co)- 
Let uy/KB{x;h) be the approximate solution given in fl8.29l) . (18.301) and let u = 
7r/i^(XMBKw) be the corresponding exact eigenf unction, where 7 = dD{fi, 2^)- 

Theorem 9.1 a) Outside any h-independent neighborhood of 0, we have 

u, BhDu = 0(e~i/(^^)) 

in L'^-norm. 

b) There exists a neighborhood Q of 0, where 

m(x; /i) = (a + r)e-'^+(")/\ (9.27) 
\\r\\mn), \\BhDr\\L2^n) = 0{h°°). 

Proof Apply flOGl) with v = u, z = X{h), e = A/i, A > 1, to get 

h\\e^^/''u\\ + h^\\BhD{e^'^''u)\\ < 0{h~^), N = N{A), 
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where we also used that e"^-!^ = 0{h-^) in ndCoe). Here, ip, = ^(0 - eg{G)) + 
0{Se) is larger than a positive constant outside any fixed neighborhood of 0, so 
u = 0{e~^/^G^^) in L^-norm there. Moreover, since ||e'^"/'^'u|| = 0{h~^), we have 

0{h-^) > h^\\BhD{e^'/^u)\\ > h^2\\e^^/^BhDu\\ - 0{h^)\\\ViJe\e^^^^u\\. 

Here Vipe = (as we shall see more in detail below), so \\e'^''^'^B^hDu\\ = 

0{h~'^~^), so B^hDu = 0{e~ch) in L^-norm away from any given fixed neighbor- 
hood of 0. The proof of a) is complete. 

To prove b), we apply fl9.26p to f = -u — X'^wkb, z = \{h) with e = Ah, A ^ 1. 
Since ||m - x^wkbH = by (EHS), and \\e^^/''{P - z)xuwkb \\ = 0{h°^) if we 

arrange so that ip^ < (J) on supp Vx, we conclude that \\e^'^'^{u — X'"WKB)||nG(Coe) = 
0{h°°) and hence 

h\\e^^/\u-xuwKB)\\+h--\\BhD{e^^/\u-xnwKB)\\ < 0{h^). (9.28) 
In a small neighborhood of 0, we have 

where by ( 19.14^ 

^-Ag{G) ^ o{l)e-^ for G{x) < C^e, 

^-Ag{G) ^ G{xY^ for G{x) > Coe. 

It then follows from (19.281) . that u = (a + r)e~^+^'^ in a neighborhood of 0, with 
MlHu) = 0{h^). 

To get the corresponding bound on B^hDr, we just have to proceed as in the 
proof of a) and use that Vipt = V0 — eg'{G)VG, where 

g\G)\/G = C(e"5) for G{x) < Got, g'{G)\/G = ^ = C>(e~^) for G{x) > Got. 

G 

Thus VV^, = V(j) + C»(v^) = 0(1) and we conclude that \\B^hDr\\L2^n) = 0{h°°). 
□ 

Remark 9.2 If we drop the assumption (19.11) and allow A^ — 1 more points p2, Pn 
in C, then Theorem 19.11 is still valid, provided that all pj^k in 18.241 with j > 2 are 
different from the value associated to pi = (0, 0). 
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10 Supersymmetric approach 



The Witten approach has been independently extended to the case of non-elhptic 
operators hke the Kramers- Fokker- Planck operator in |22j (in supersymmetric lan- 
guage) and in |2j (in terms of differential forms). See also [18] . 

We start by a quick review of that in the semiclassical case, then we establish 
some basic facts about the principal and subprincipal symbols, especially at the 
critical points of the given weight function. 



10.1 Generalities 

Let 

A{x) : T*M T^M, x e M, (10.1) 

be an invertible map depending smoothly on x G M. Then we have the real non- 
degenerate bilinear form 

{u\v)a(x) = {A^A{x){u)\v), u,ve A^T*M. (10.2) 

If a : A^T*M A^T*M is a linear map, we define the "adjoint" a^'* : A^T^M 
A^T*M by 

{au\v)A(x) = {u\a'^'*v)A{x) (10.3) 

(In the complexified case, we use the sesquilinear scalar product {u\v)a = {u\v)a 
and define a^'* the same way.) 

If is a one form and u and v are k — 1 and k forms respectively, we get at x: 

{uj A u\v)a = {Auj A {a''-^A)u\v) = {A^-^Au\{Auyv) = {u\{Aujyv)A, (10.4) 

so 

(cu^)^'* = (Acj)J, (10.5) 

where J denotes the usual operator of contraction. Let n{dx) be a locally finite 
measure with a smooth positive density. When M = R", /i will be the Lebesgue 
measure. Sometimes we also use the symbol /i for the corresponding density. 
If u,v are smooth k forms with suppw fl suppf compact, we define 

{u\v)a = j {u{x)\v{x))A{x)^'{dx), {u\v)a = j {u{x)\v{x))A{x)^^{dx) 

and denote by a^'* the formal adjoint of an operator a : C^{M]A^T*M) 
P'(M; A^T*M). If we fix some local coordinates xi, and write fi{dx) = ^{x)dx 
(by slight abuse of notation), we can consider 

dx^ : C^iM; A^T*M) C^{M] A^T*M), 
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acting coefficient- wise, and a straightforward computation shows that 

(5.,)^'* = {fiXA'A))-\-d.,^) o ifiXA'A)) = -d,^ - ^ - tA'=A)-iH^,(AM)). 

(10.6) 

We only retain that 

ihd,^)^'* = -hd,^ + Oih), (10.7) 

where 0{h) stands for multiphcation by a smooth matrix, which is 0{h) with all 
its derivatives, uniformly on all of M when M = R", and which is = when A{x), 
n{x) are constant. 

Let (j) e C°°{M] R) and introduce the Witten (de Rham) complex 

= e-'^/^ohdoe'^''' = hd+{d(t))'' : C^{M- a''T*M) C^{M; a''+^T*M), (10.8) 

with ^0 = 0. 

In local coordinates (always the canonical ones when M = R") we have 

n 

d^ = ^(/la,,. + d,^<P) o dx^, (10.9) 
1 

where hd^^ + c?xj</' acts coefficient-wise and commutes with dx'j, so 

n 

<'* = J2^-hd., + d^J + 0{h)) o Aidx^y , (10.10) 
1 

where from now on, 0{h) and 0{h'^) will have the same meaning as after fll0.7p . 
The corresponding Witten-Hodge Laplacian is given by 

-AA = df*d^ + d^df*. (10.11) 

SincG 

we also have 

-AA = {d^ + dfy, 

and —Aa conserves the degree of differential forms. 

Choose local coordinates xi, ...,Xn (to be the standard ones when M = R") and 
write 

A{dxk) = ^Aj^kix)d^^, Aj^k{x) = {A{dxk)\dxj). 
j 
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Let Zj = hd^. + d^.cj), Z^'* = hd^^* + Notice that [Zj, dx^] = 0, so 

[Zt'\Aidx,y] = 0. 
Writing = Xli Zj o dx^, df* = X]" Z^'* ° M^Xk)^ , we get 

-A^ = J2{Z^'* o Aidxk^dx^Zj + Zjdx^AidxkYz^'*). 

Here, we use the general identity u'^jj,^ + jjL^v^ = {u, n)l on the first term in the 
parenthesis to get 

-Aa = I + 11 + III, 
I = '^Z^'*Aj^kZj, 

II = -J2zt'*dx^MdxkyZj 
111 = ^dx'^ZjZ^'*A{dxk)K 
where Aj^k = {dxj, A{dxk)) ■ We have 

[z„A{dxky] = h{{d^^A){dxk)y. 

Using the identity (f/^.*)'^.* = JJ (see Subsection I10.4p . we see that 

[Zt'*,dx^t^'* = [CA{dx,)y,Z,] = -hi{d,^'A){dx,)y, 

hence 



[Zt'*,dx^] = -hCA~'d^Mdxj)r 



Using these commutator relations, we move dx^ to the left and A{dxk)^ to the right 
in II and combining with III, we finally obtain that 

-A^ = J2^k'*A,,,Z, + J2dx^[Z,,zt'*]Aidx,y (10.12) 

j,k j,k 

+hJ2zt'*dx^{id,^A)idxk)y + hJ2C^~\d.jA)dx,)''A{dxkyZ, 
+/i2^(*A-\9,,*A)cix,)^((a,^A)dxfc)J. 

j,k 



59 



Modulo 0{h){hd^ + d^cj)) + 0{h){-hd^ + d^cfy + ^(/i^), we get 

-A^ = Y^{-hd,^ + d^^<P)A^,k{x){hd,^+d^^(P) (10.13) 

where the error terms vanish when A{x) and are constant (for the chosen 

coordinates). 
Now write 

A{x) = B{x) + C{x), 'B{x) = B{x), tJ{x) = -C{x). (10.14) 
Then f lTIU3D gives 

-A^ = J2hhd.,, + d,J)Bj4x){hd,. + d,^^) (10.15) 

j,k 

+ ^ii9xk(l>)Cj,khd^^ + hd^^ o Cj^k o (5xfe0)) - ^ hdx^{Cj^k)hdx^ 

j,k j,k 

+j2'^hd,^d,jodx^A{dxky. 

Again, this becomes an equahty when A, ^ are constant. Note that the last term 
vanishes on 0-forms, i.e. on scalar functions. To recover the Kramers-Fokker-Planck 
operator (cf [22]), replace n by 2n, put M = Ri^" 



2 \-l 7 

and let fi = dxdy be the Lebesgue measure. Then fllO.lSp is an equality and we get 
for 0-forms: 

n 

-A^°) = lY^{-hdy^ + dy^mdy, + dy^<j>) + (IO.I6) 

+- ^[idy,(l))hdx^ - {dx^(f)) hdy^ + hd^^ o - hdy^ o d^^cf)] 
j 

n 
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where 

is the Hamihon field of 4> with respect to the standard symplectic form ^ dyj Adxj. 
If we choose 

<P{x,y) = ^i/ + V{x), (10.17) 
we get the Kramers-Fokker-Planck operator 

- Af =yhd,- V'{x) ■ hdy + |(-H + y) ■ (hdy + y). (10.18) 

10.2 The principal symbol of the Hodge Laplacian 

The principal symbol of — in the sense of /i-differential operators is scalar and 
given by 

j,k j,k 

The corresponding real symbol q{x,^) — —p{x,i^) is given by 

q{x,o = 5^^,-fe(a + a,,0)(e,-a,^.0) (10.20) 
j,k 

j,k j,k 

It vanishes on the two Lagrangian manifolds A±^. 
We define 

i^± = H . (10.21) 
Using xi, Xn as coordinates on A±<^, we get 

I/+ = 2j2Aj,kdx,<pdx, = 2A{x){(l)'{x)) ■ dx (10.22) 

j,k 

u_ = -2Y,Aj,kdx.<pdx, = -2'A{x){<f>'{x)) ■ dx. (10.23) 

j,k 



61 



(Even more radically, we could say that = 2A{x){d(j){x)), where A{x) is viewed 
as a map T*M — >■ T^M, and similarly for z/_.) 

Let Xq be a non-degenerate critical point of 0, so that and A_0 intersect 
transversally at (a;o,0). The spectrum of the linearization Fg of Hg at (a;o,0) is 
equal to the union of the spectra of the linearizations 

z/° = (2A(xo)0"(xo)x) ■ and z/° = -(2*A(xo)0"(xo)x) ■ (10.24) 

of and z/_ respectively at Xq. Thus we are interested in the eigenvalues of the 
matrices A(f)", ^A(f)", where we write A = A{xo), and 0" = (f)"{xo) for short. Here, 
we notice that = (j)"~^\A(j)")(j)" has the same eigenvalues as Acj)" and similarly 
(p" A, (j)"*A are isospectral to A(j)". Thus 

The eigenvalues of Fg are given by ± 2Aj, (10.25) 
where Ai, A„ are the eigenvalues of Acp" . 

From questions about hypoellipticity (see [23]) we would like to know when all 
the eigenvalues of Fp avoid the real axis, or equivalently, when all the eigenvalues of 
Fg (the linearization of Hg at (xo,0)) avoid the imaginary axis. 

We assume from now on that 

B{x) > 0, X EM. (10.26) 

Then, if 0°(a;) = ^(f)"x ■ x is the Hessian quadratic form of at xq, we have 

zyO(0O) = {2B(j)"x,(f)"x) > 0. (10.27) 

Lemma 10.1 Let fi{x, dx) = Mx-d^ be a real linear vector field on R". Let n± G N, 
+ n_ = n. Then the following two statements are equivalent: 

(A) M has n+ eigenvalues with real part > and n„ eigenvalues with real part < 0. 

(B) There exists a quadratic form G : R" — > R of signature {n^,n_) and a constant 
C > 0, such that 

fi{x,dx){G) > X e R'^. (10.28) 

Proof Assume first that (n+,n_) = (n,0). If (A) holds, we know that 

||e**^xf > ie*/^||xf , t > 0, 
for some constant C > 0, and we can put (by a classical argument) 

G{x) = Grix) = ^ f We'^^xfdt, T > 1. 
T Jo 
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Then 



Mx ■ d^iGrix)) = ^ [ Mx- d^\\e'''xfdt = ^ / 



dr 



= ^(l|e™xf - llxf ) > i(ie^/^ - l)||xf > ^||xf , 

if T is large enough. Thus we get (B). 

Conversely, if (B) holds, we have with a new constant C > 0, that 

n{x,da:)G > ^G, 

and hence 

^G{e'''x) > ^G{e'^'x), 
dt C 

so 

G{e'^^x) > e*/^G(a;), t > 0. 

Thus ||e*^-^x|| > C-ie*/^||x|| for some new positive constant C and we conclude that 
the eigenvalues of M all have positive real parts. 

Now consider the case of general (n+,n_) and assume first that (A) holds. Then 
we have the M- invariant decomposition, R" = © where dim {L±) = n± and 
a{M\^^) belongs to the open right half plane in the + case and to the open left half 
plane in the — case. Hence we have positive definite quadratic forms G± on L± such 
that 

±Mx ■ > X e L±. 

Then G = G+ © {~G^) (defined in the obvious way) has the required properties in 
(B). 

Conversely, assume that (B) holds. Let L+ be an n+-dimensional subspace on 
which G is positive definite. By (110. 28p . we have for all x G R", 

^I{x,^,)G{x) > ^G{x), OO, 

so if X G we get 

G{e'^'x) > e*/^G(x), t > 0, 
and hence with a new constant C > 0, 

|e**^a;| > ^e*/^|x|, t > 0. 
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Similarly, if dim (L_) = n_ and G is negative definite on L_, we get 

|e**^x| > ^el*l/^|x|, X e Z_, t < 0. 

Now we have the M-invariant decomposition 

R" = L+ ©Lo©L_, 

where L^, L^, are the sums of generalized eigenspaces of M corresponding to the 
eigenvalues with real parts > 0, = and < respectively. We see that necessarily, 
L+ n (Lo © L_) = 0, so 

n+ = dimL_(. < dimL+. 

Similarly, 

n_ = dimL_ < dimL_, 
so dimL+ = n+, dimL_ = n^, dimLo = 0, and (A) follows. □ 

If 5 > 0, then (110.271) with strict inequality for x ^ 0, and Lemma [10.11 imply 
that z/^ has n± eigenvalues with ± real part > 0, where (^4., «_) is the signature of 
(j)"{xo). It follows in that case that Fp has no real eigenvalues. This last conclusion 
also follows from [23]. Indeed, in that case the quadratic approximation p^{x,C,) of 
p at (xo,0) is elliptic in the sense that \p^{x,^)\ ^ |xp + and takes its values in 
an angle — | + e < argp° < | — e for some e > 0. 

Now return to the general case, when we only assume (110.261) and 0" = 0"(xo) 
is non-degenerate of signature (n+,n_). 

We next make some remarks about the quadratic approximation of p at (a;o, 0). 

Proposition 10.2 a) Assume that the matrix A(f)" of z/° has m± eigenvalues with 
dz real part > 0, m+ + m_ = n. Then there exists a real quadratic form Q{x,C,) on 
R^" such that 

Rep°((x,0+^effg(x,0)>^|(a:,OI', (a:,OeR'", < e « 1. (10.29) 

b) Conversely, assume that there exists a quadratic form Q such that (110.290 holds. 
Then Acj)" has n± eigenvalues with ± real part > 0, where {n+,n-) is the signature 

of no). 
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Recall that the condition fl4.2ip implies the existence of Q as in a) of the propo- 
sition. (The converse is not true however. It is easy to find examples of purely 
imaginary quadratic forms for which there exist Q as in a) of the proposition.) 
Proof a) Choose of the form G{x) + so that Hg{x,i) = (G'^, -G'J. 

Recall that p° is the quadratic approximation of p at (xq, 0), obtained from (110. 19p 
by freezing ^ at Xq and replacing dx(p by (j)"x = (0°)'(x), 0" = (f)"{xo). 

We get 

Rep°(x,0 = $^i?,,fc(e,e;^ + 5.,0°5.,0°) 

+2e^Y.C,Ad.J>.,G-Y,C,A<P''d^GU,)+Oie\x,0\^^^ 

j,k j,k 

SO it suffices to have 

C(j)"x-d^G > hx\^ onf-'XiB), (10.30) 

cj)"C^-d^G > h^\^on^^{B), (10.31) 

where we also used the antisymmetry of C in the last equation. This would follow 
from 

A(l>"x ■ d^G > ^\x\\ ■ d^G > ^leP, G R^ 

and in order to find such functions G and G, it suffices to apply Lemma [10.11 to the 
isospectral matrices A(f)" and (j)"A. 

b) Let Q be as in (110. 29p . The quadratic form 

p'{x, = /(exp teHgix, 0) = /(x, + teHgp' + ^(e^) 

is elliptic on the real phase space and takes its values in an angle 

■\ 7T < e TV e 1 r r 

e*l-2 + c'2-cl[0,+oo[, 

so we know from [23] that FpO has no real eigenvalues. 

On the other hand p^ and p° are related by a canonical transformation, so FpO 
and Fpg are isospectral. Hence A(f)" has m± eigenvalues with ± real part > 0, where 
m+ + m_ = n. To see that m± = n±, we just replace -B by -B + 51, < 5 <^ 1, 
to reduce ourselves to the elliptic case, and apply the observation after the proof of 
Lemma [TIUl □ 
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10.3 The subprincipal symbol 



We next look at the subprincipal term in (110.131) . Write 

A{dxk) = ^Ay^kdx,, A(rfxfc)J = ^Ay^kdi^, 

so the second sum in (110.131) becomes 

j,k,v j,u 

which simplifies further to 

2hJ2i<P"°'A)idx,)%. (10.32) 

j 

Now we restrict the attention to a non-degenerate critical point xq of (p and 
we shall compute the subprincipal symbol of —Aa at the corresponding doubly 
characteristic point (xo,0). At that point (p" o */l : T*^M —>■ T*^M is invariantly 
defined and it is easy to check that (110.321) is also invariantly defined: we get the 
same quantity if we replace dxi, ..,dxn, d^^, .■,dri;„, by uji,..,ujn, ^i,--,^ny where 
LJi, ..,Un is any basis in the complexified cotangent space and ujI, ..,u!^ is the dual 
basis of tangent vectors for the natural bilinear pairing. 

Assume that the equivalent conditions of Proposition 110.21 hold and denote the 
corresponding eigenvalues (that are also the eigenvalues of 0" o *A) by Ai, .., A„ with 
Re Xj > for 1 < j < and with Re Xj < for n+ + 1 < j < n = + n_ . The 
eigenvalues of Fp are then ±2iXj (in view of (110. 25p and the isospectrality of Fp and 
iFg reviewed prior to Lemma IHTTl) , so 

^ n+ n 

^^Pp--=- E /^ = E2^^- E2\- (10.33) 

Im /i>0 

The subprincipal symbol of the first term in (110.131) (at (xq, 0)) is equal to 
1 " 

j,k j,k 1 

(10.34) 

The eigenvalues of ^j{(t>" ° ^A){dxj)^dij on the space of m-forms are easily cal- 
culated, if we replace dxi, .., dxn by a basis of eigenvectors a;i, ..,a;„ of so that 
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and by the corresponding dual basis vectors u*. (Here we assume to start 
with that there are no Jordan blocks. This can be achieved by an arbitrarily small 
perturbation of A, and we can extend the end result of our calculation to the general 
case by continuity.) We get 

3 3 

A basis of eigenforms of this operator is given by uJj-^ A .. A Uj^, 1 < Ji < J2 < •• < 
jm < n and the corresponding eigenvalues are Aj^ + .. + Xj^. 

Let Sp be the subprincipal symbol of (110. 13p at (xq, 0). Then the eigenvalues of 

-trFp + Sp, acting on m forms 

are 

n+ n n 

E - E - E + 2(^.1 + •• + V) (10-36) 

1 n++l 1 

n 

= 2(Aj, + .. + Xj^ - E 1 < Ji < •• < Jm < n. 

n++l 

We conclude that if m 7^ n_, then all the eigenvalues have a real part > and if 
m = n_, then precisely one eigenvalue is equal to 0, while the others have positive 
real part. 

10.4 A symmetry for adjoints 

Our last remark in this section concerns symmetry relations for the A, * adjoints. If 
D : L^{n; A''T*n) L'^{n; A^T*n), then a simple calculation shows that 

D^'* = (a'^(U))-^D*(A^'(U)), 

where D* denotes the adjoint with respect to the measure /i. 
We also have 

{u\v)a = {v\u)tA. 
Playing with these relations we see that 

D = {D^^*y^'*. 

This can be applied to — and we get 
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11 The double well case 



In this section we assume that 

M = R", A = A{x) is independent of x, and invertible. (H-l) 

We decompose A as in (fTOlD and assume (^M^. Let G C°°(R";R) be such 
that 

= ((59,0,9,0)) = 0(1), |a|>2. (11.2) 

Consider P'^'^^ = — A^'' which according to fllO.lSp becomes 

p(i) = J2hD..B,,khD,^ + J2i9x,(i))BjAd.J)-htT{B<j)") (11.3) 

j,k j,k 

+ J2^{d.J)C,,khd,^ + hd,^ o C.Ad.J)) + 2h 5^(9,,9,,0)rfa;;A(rfx,)J. 

j,k j,k 

Apart from the third and the last terms which are 0{h) with all their derivatives, 
this is of the form (12.11) with bj^k = -Bj,fc, Cj{x) = J2k Cj,kdxk(p, Po{x) = {Bd^cp, <9,0). 
We define P2,Pi,Po as in Section [2] and see that fl2.4l) - fl2.7p hold. 
Assume that 

is a Morse function with critical points xi, a; at G R", (11-4) 

|0'(x)| > 1/C, \x\ > C. (11.5) 

Let C = {pj] j = 1,...,A^} where pj = {xj,0) G T*R". Then pj are critical points 
with critical value for p2,pi,po- Since A is invertible, Mi^B) nA/'(C) = 0, and hence 
(14.41) holds. We adopt the dynamical assumptions (14.211) . (I4.22p (or equivalently 
(5211), (1121) and (031) ). Then we can apply the results of Sections El E] to p('') 
since neither the presence of the bounded subprincipal symbol in (I11.3P nor the non- 
scalar nature of the operators has any serious influence. In the preceding section 
we saw that we are in the case when the conditions of Proposition 110.21 are fulfilled. 
The only change in Theorem 18.31 is that the pj^k in (18.241) are of the form 

1 " 1~ 
= - X] ^j.kA^ + {Fp{pj)) + 7i,fc, (11-6) 

k=l 

where ■jj^k is any eigenvalue of the subprincipal symbol Sp{q) at (xj,0). From the 
calculations in Subsection 110.31 we notice that the pj^k will be confined to a sector 
{0} U {|argz| < tt/2 — 1/C} around [0, +oo[ and it is precisely when Xj is of index 
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q (i.e. when the Hessian of at Xj has precisely q negative eigenvalues) that one of 
the fij^k may be equal to 0. 

We now add more specific conditions for the double well case. Assume that 

(p has precisely three critical points, two local (H-''') 
minima U±i, and a "saddle point" Uo of index one. 

Then (f){x) +00 with 0(a;) > ^\x\, for |x| > C. 

Put Sj = (j){Uo) - (f){Uj), j = ±1, so that Sj > 0. The set (f)^\] - oo,0(f/o)[) 
has precisely two connected components Dj, j = ±1, determined by the condition 
Uj G Dj. Under these assumptions we know that P'^^-' = — has precisely two 
eigenvalues /xq,/^! = o{h) spanning a corresponding 2-dimensional spectral subspace 
E^^\ Actually one of these two eigenvalues, say fxo, is equal to with e"'^/'' as the 
corresponding eigenfunction and since a truncation of this function can be used as a 
quasimode near each of U±i we also know that /ii = 0{h°°) (cf (111.131) ). Moreover, 
— A^'* has precisely one eigenvalue Jli = o{h) and — A^^ has no eigenvalues = o{h) 
for k > 2. Since our operators are real we know that the spectra are symmetric 
around the real axis, hence /io,/^i,/^i are real. From the intertwining relations 

-A^)rf, = d,i-Af), -Aj)<- = -<'*A«, 

we then also know that Jli = fii. 

In fact, when i? > it follows from the ellipticity and the estimates in Section 
E] that all eigenforms and generalized eigenforms corresponding to an eigenvalue in 
D{0,Ch) belong to 5(R"), so if fii ^ and {-Af - /ii)u = 0, u G L^, then 
u e iS(R") and 7^ d^u G 5 is an eigenform for — A^"* with the same eigenvalue. 
A priori we cannot exclude that /ii = and that —A^^^u = Const e"'^/''. Then 
again, 7^ d^u G 5 is a corresponding eigenvector of — A^"*. In the general case, 
we let < ^ B =: Bq when e \ 0. (Take for instance B^ = B + el.) On a 
small circle -D(0, h/C) we know that {Pe'^^ — z)^^ = C(^) uniformly for < e ^ 1. 
If {Pq'^^ — z)u = V, u,v G S (for e = 0), then {Pj:'^^ — z)u = v + r^, ^ 0, so 
(Pi"^ - zy^v = u- (Pi''^ - z)-V, ^ n, e ^ 0. Since (P(^) - z){S) is dense in we 
conclude that (Pi"^"* — z)^-^ (P*^'^^ — z)^^ strongly for \z\ = h/C. We have then the 
corresponding fact for the finite rank spectral projections and their compositions 
with Pi''^ 

— / {z- P.y^dz and — / z{z- P,)-^dz. 
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These are finite rank operators and converge therefore in norm. It follows that 

a(Pi'')) n D{0, h/C) ^ (T(P('?)) n D(0, h/C), when e ^ 0, (11.8) 

so we get /ii = jii also in the general case. 

Let Xj £ C^{Dj) be equal to 1 on Dj fl — oo, (j){Uo) — eo]) for eo > fixed 
but arbitrarily small. Consider 

/, = /i-"/%(/.)e-^W^')-*(^^))x.(a:), J = ±1, (11.9) 
where cj ~ Cj^ + /iCj^i + ... > is a normalization constant with cj^ > 0, such that 

||/,|| = 1. (11.10) 

We also have 

p(o)(/^.) = [p(o),;^^.](c,-e-^('^(")-<^(^^«) = 0(/i-^Oe-^(^^-^o)), (11.11) 

for some A'^q > 0. 
If 

n(°) = ^ ^(^ - 7 = 9^(0, ^) (11.12) 

is the spectral projection of P^^^ onto i^'-^-' we know from Theorem 18.41 that U^^^ = 
0{1). It follows from ffTLTTll that 

e, := n(°)/i = fj + Oih-^'e-^^^^-''^) in L\ (11.13) 
In fact, we write (111. Ill) as P^'^^ fj = rj, 

{z~P^''>){f,) = zf,-r^, 

{z - P(°))~V, = -/,• + {z- P^^'^Y^z-^j 

and integrate, using the bounds on the resolvent provided by Theorem 18. 4[ 
From (111.131) we see that 

||e,f = 1 + C(/i-^2e-^(^^-^°)), (11.14) 

(ei|e_i) = 0(/i-^^e-^(^-'°-^o)), (11.15) 

where 

Smin = min(S'„i, ^i). (11.16) 
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Let be the one-dimensional eigenspace of P^^-* corresponding to /ii. From 
an easy extension of Theorem 19.11 to the non-scalar case with the presence of other 
non-resonant wells {U±i) as in Remark 19.21 we know that E^^^ is generated by an 
eigenform 

eo(x; h) = xo(a;)e-^'^+(^)/i-?ao(x; h) + C(e-^«/''), (11.17) 
where Xo ^ C^(neigh (Uq)) is equal to one near Uq, Sq > 0, 

oo 

ao{x; h) ~ y^^aQ,k{x)h^ 



is a symbol as in Theorem 19.11 with ao,o(f^) 7^ 0, and 0+ G C°° (neigh (?7o); [0, oo[) 
satisfies 

0+(x) - |x- [/o|^ (11.18) 

and solves the eikonal equation 

g(x,0'+(x))=O, (11.19) 

with q = p2+pi -po, 

P2 = {Bix)^,0, Pi{x,0 = 2{Cix)<P\x),0, Poix) = (5(x)0'(x),0'(x)). 

A^^ is the stable outgoing manifold through {Uq, 0) for the Hg-fiow and recall 
that 0'^(t/o) > by Proposition 18. 2[ (Similarly we have a stable incoming manifold 
A(^_.) Let k± be the number of eigenvalues of the linearization of Hg^^ at that point 

with ± real part > 0, so that /c+ + = n. Let K^, C be the corresponding 
stable outgoing and incoming submanifolds of dimension k+ and k- respectively. 
Then C A^^, C and — 0(?7o) — 0± vanishes to the second order on 
T^x{.K±). Since (p"{Uo) has signature (n — 1, 1), we conclude that dimi^^ = n — 1, 
dimK_ = 1. (This also follows from Proposition 110.21 ) It is also clear that A,^, A,^^ 
intersect cleanly along K±, so we get 

0+ - (0 - 0(f/o)) - dist (x, 7r,.(ir+))2, - 0(f/o) - 0_ ^ dist (x, 7r,.(ir_))2. (11.20) 

We next make some remarks about the adjoint operator — At4 = (— A^)*^'* (cf 
Subsection 110.41) . The principal symbol is p2 — ipi +Po = p{x, —$,) = p{x,^) and the 
corresponding real "g"-symbol is q{x,^) = q{x, — ^). Since our dynamical conditions 
are invariant under a change of sign of the ifp^ -direction, all our assumptions are 
equally valid for — At4. This also holds for the geometric discussion above, so if A,^;^^, 
A^*_ denote the outgoing and incoming if^-invariant Lagrangian manifolds through 
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(f/o, 0) and K"^ C A,^ the outgoing/incoming manifolds for Hq, (noting that g = 
on A(^), then dimi^^ = — 1, dimi^'l = 1 and 

0; - (0 - 0(f/o)) - dist (a;, v^,(K:))^ - 0(f/o) - 0*_ - dist (x, vr,(ir:))2. (11.21) 

In view of the general relation 

J,{Hq) = -Hg, where J : (a;, ^ (x, -0, (11-22) 

we see that A^;^ = J(A0^), A^;^ = J(A^_), or more simply 

0*_ = -0+, 0: = -0-, (11.23) 

giving in particular from flll.2Up . flll.2ip . 

- 0(?7o) + 0; - dist (x, 7i,{K_)f, - 0(f/o) + 0+ - dist {x, 'K^{K*_)f. (11.24) 

Let /Iq = 0, /i^ be the two eigenvalues of P^^^ := — Aj^-* that are o(/i) and let 
ni°^ be the spectral projection onto the corresponding spectral subspace Ef\ Then 
e* = nl°^/j, j' = ±1, span and satisfy (111.131) . Similarly for P^^^ = — At^-* we 
have the generating eigenform 

e*(x; h) = xo(a;)e-^'^+(")/i-?a*(x; h) + C(e-^«/'^) in (11.25) 

for the one dimensional eigenspace Ei^^ corresponding to fi^. 

Now, using that our eigenvalues and operators are real, we know by duality that 

/it = /ii, (11.26) 

and that (^f^E^) and {Ei^\ E^^^) are dual pairs for the scalar products {u\v)l^ 
and {u\v)a respectively. In fact, {{z - Pi^V^)^'* = {z - P^''^)-\ 

From Subsection 110.31 we know that ao,o(t^o) is an eigenvector corresponding to 
the negative eigenvalue of 0"o*A at Uq, and aQo(t^) is an eigenvector corresponding 
to the negative eigenvalue of 0" o A. Since (0" o A)^'* = 0" o *A, we know that the 
two eigenvalues are equal and that the A-product of the two eigenvectors is 7^ 0; 

(a*^o(f/o)|ao,o(f/o))A^O. (11.27) 

It follows that (eQ|eo)A ^ 1 and after renormalization of Cq we may assume that 

{e*o\eo)A = 1. (11.28) 
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Similarly, using (111.101) 

{e]\eu)=5,^k + 0{e-^), j, = ±1. (11.29) 

Let ( A_i Ai ) be the matrix of : E^^^ — >• E'^^'> with respect to the bases 
(e_i, ei ) and (cq). (Strictly speaking, we approximate our operators by elliptic 
ones as in flll.Sp and pass to the limit.) Let 

K 

be the matrix of d^'* for the same bases. The eigenvalue /ii can be viewed as the 
second eigenvalue of d^'*d^ : E''^^ — > E^^^ or equivalently as the scalar d^d^'* : 
(using also that P*^^^ has no eigenvalue = o{h)). Either way, we get 



We get 



where 



^1 = AliA_i + A*iAi. (11.30) 

Afe = {e*Q\d^ek)A, k = ±1, (11.31) 

X* = {9j\df*eo)A, J = ±l, (11.32) 

9-1 9i) = {eU el){l + 0{e-^)) (11.33) 



is the base in Ei^^ that is dual to (e_i Ci ). Here the complex conjugate signs are 
superfluous since we work with real operators, eigenvalues and functions. 

Let X £ C^(neigh (Uq); [0, 1]) be equal to 1 near Uq. Using that (i^C-i = A_ieo, 
we get, dropping the bars from now on, 

A_i = (e*|V-i)A (11.34) 
= ie*o\xd<pe-i)A + X-iiel\{l- x)eo)A 
= (e*|xc/0e_i)A + C(e-^)A_i. 



Here 



iel\xd^e^i)A = ie*Q\[x,d^]e^i)A + {d'^'*e*o\xe-i)A- (11.35) 



Now the matrix of d^'* : Ei^^ Ei^^ with respect to the dual bases is the adjoint of 
the one of d^ : E^^'^ E''^\ so d^'*eQ = X-19-1 + Xi9i, and expressing gj as linear 
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combinations of the e* by means of flll.33p and using (111.131) for the e^^ we see that 
the last term in flll.24p is of the form 

idf*el\xe^^)A = 0(e-c^)A_i + 0{e'^)\,. (11.36) 

Thus we have obtained 

(1 + 0(e- ck))A„i + 0(e-c^)Ai = (e*| [x, d^]e.,)A = -h{el\{dx) A e.i)^, (11.37) 

and we shall study the last expression. The contribution from the remainder in 
(nrm is 0{h-^^)exp}^{-S-i + eo - ^) = C(l)e-^(^-i+2Z') if we choose eo small 
enough. A similar estimate holds for the contribution from the remainder term in 
(lll.25p . As we shall see, the contribution from the leading terms in f lll.l3p . flll.25p 
will be larger. It is equal to 

-c.i{h)h'-^ J X-i(a;)(A(x)a*(x;/i)|rfx(2;))e-^('^+(")+*(")-'^(^-^»(ia:. (11.38) 

Here by ^LM, 

0;(x)+0(x)-0([/_i) =0;+(0(x)-0(t/o))+5_i ^5_i+dist(x,7r,(K_))2, (11.39) 

so we expect flll.38p to behave like some power of h times e~T^^-'^ with the main 
contribution coming from a neighborhood of supp (dx) H supp (x-i) H nx{K^). Now 
0(x) — 4>{Uo) ^ —\x — f/oP on iTx{K^) while x~i has its support in D_i and equals 
1 in the subset of D_i where 0(a:) — 0(t/o) ^ with eo arbitrarily small. Because 
of the presence of dx which has its support in an annular region around Uq, we see 
that X-i{^) = 1 in neigh (-D_i fl supp (dx) H 7r^(-ft'_)), so we can forget about X-i in 
(I11.38P and just integrate over neigh (supp (dx) H Hx{K^i), -D-i). 

Let us look at {A{x)alQ{x; h)\dx{xi)) at a point xi G supp (dx) nnx{K^). At Uq 
we know that a^ Q is an eigenvector of (j)" o A associated to the negative eigenvalue, 
so Aal Q is a corresponding eigenvector for Acf)" which is the linearization of ^Hg, 

Hence Aa^Q at Uq is tangent to T^^iK-). This will remain approximately true at 
Xi since the latter point is close to Uq. Choosing x to be a "circular" standard 
cut-off, we see that {AaQQ,dx) is non-vanishing of constant sign at every point 
Xi G 7r^(fC_) n D^i where dx 7^ 0. 

By stationary phase it is now clear that the integral (111.381) is equal to 

/i^£„i(/i)e-T^^-i, £_i ~£_i,o + /i^-i,i + ..., ^-1,0 7^0. (11.40) 
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Returning to (I11.37P and modifying i_i by an exponentially small term, we get 

(1 + C(e-^))A_i + C(e-^)Ai = hk^i{h)e-^^-\ (11.41) 

Similarly, 

0(e-^)A_i + (1 + 0(e-^))Ai = hkiih)e-^^\ £i ~ ^i,o + + ii,o ^ 0. 

(11.42) 

Inverting the system, we get 

(-.)Hi.o(e-^))('tf-;{tl::j::). (u,43) 

Now turn to A* in flll.32p . In view of flll.33p . we have 

(^;.')^(l + 0(e-A))(-.). (U.44) 

where 

aj = {d^e*\eo)A = {eo\d^e*)tA 

which can be identified with the expression (111.310 after replacing A by M and 
making the corresponding substitutions, Cq — > Cq, Cj —>■ e*. Hence we have the 
analogue of ( 111.43^ . 

e*{h)r.ei, + hq, + ..., £*o^o. 

We finally claim that ^j,o^jo > 0- Indeed, this number is real and different form 
zero and if we deform our matrices to reach the self adjoint case (with A > 0) we see 
that we have a positive sign). 

Combining this with (lll.30p we get the main result of this work: 

Theorem 11.1 Let P = -A^^^ where we assume (fTLTj) . (nT2^ . ffTLTD . We also 
assume that P satisfies the additional dynamical conditions fl4.2ip . (I4.22p (or equiv- 
alently fl4.2ip . (I4.24p ) and (14.230 . Then for C > large enough, P has precisely 2 
eigenvalues, and fii in the disc D{0,h/C) when h > is small enough. Here fii 
is real and of the form 

/ii = h{ai{h)e-^^'/'' + a_i{h)e-^^-'/^), (11.46) 

where aj{h) are real, aj{h) ~ ajfl + ajih + a^- o > 0, Sj = (f>{Uo) — (piJJj). 
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